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Abstract. Particle number fluctuations and correlations in nucleus- nucleus collisions 
at SPS and RHIC energies are studied within microscopic transport approaches. In 
this review we focus on the Hadron-String-Dynamics (HSD) and Ultra-relativistic- 
Quantum-Molecular-Dynamics (UrQMD) models The obtained results are compared 
with the available experimental data as well as with the statistical models and the 
model of independent sources. In particular the role of the experimental centrality 
selection and acceptance is discussed in detail for a variety of experimental fluctuations 
and correlation observables with the aim to extract information on the critical point 
in the (T, fis) plane of strongly interacting matter. 
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1. Introduction 

1.1. Relativistic Heavy-Ion Collisions 

Hot and dense nuclear matter can be generated in the laboratory in a wide range of 
temperatures and densities by colliding atomic nuclei at high energies. In the collision 
zone the matter is heated and compressed for a very short period of time. If the energy 
pumped into the formed fireball is sufficiently large the quark-gluon substructure of 
nucleons comes into play. At moderate temperatures, nucleons are excited to short- 
lived states (baryonic resonances) which decay by the emission of mesons. At higher 
temperatures, also baryon-antibaryon pairs are created. This mixture of baryons, 
antibaryons and mesons, all strongly interacting particles, is generally called hadronic 
matter, or baryonic matter if baryons prevail. At even higher temperatures or densities 
the hadrons melt and the constituents, the quarks and gluons, form a new phase, the 
Quark-Gluon Plasma (QGP). 

High-energy heavy-ion collision experiments provide the unique possibility to create 
and investigate these extreme states of matter. The study of matter at extremely high 
temperature and baryon density, where the hadronic matter transforms to a QGP is the 
aim of a variety of experiments at current and future facilities: NA38, NA49, NA50, 
NA60 and NA61/SHINE at the Super-Proton- Synchrotron (SPS) |lH5]; PHENIX, 
STAR, PHOBOS and BRAHMS at the Relativistic-Heavy-Ion-Collider (RHIC) [BHS]; 
ALICE at the Large-Hadron-Collider (LHC) [lO]; CBM and PANDA at the Facility for 
Antiproton and Ion Research (FAIR) pTj; MPD at the Nuclotron-based Ion Collider 
fAcility (NICA) [I2]. 

Relativistic nucleus-nucleus (A+A) collisions have been studied so far at beam 
energies from 0.1 to 2 AGeV at the SIS (Schwerlonen-Synchrotron), from 2 to 11.6 
AGeV at the AGS (Alternating Gradient Synchrotron) and from 20 to 160 AGeV at the 
SPS [inilTl]. While part of these programs are closed now, the heavy-ion research has 
been extended at RHIC with Au+Au collisions at cm. nucleon pair energies -^/snn from 
~ 20 to 200 GeV (equivalent energies in a fixed target experiment: 0.2 to 21.3 ATeV). 
In the near future, further insight into the physics of matter at even more extreme 
conditions will be gained at the LHC, which will reach center-of-mass energies of the 
few TeV scale. Apart from LHC, the SPS successor SHINE operates at CERN in order 
to scan the 10A-158AGeV energy range with light and intermediate mass nuclei 0. 
At FAIR, which is expected to start operation in 2015, collisions of gold nuclei from 5 
AGeV up to 35 AGeV will be studied. At NICA it is planned to start the experimental 
program of colliding Au and/or U ions as well as polarized light nuclei at energies up 
to of 5 AGeV in 2013 (an upgrade to ^/sNN = 9 GeV is foreseen [15]). 

At very high beam energies - as available at RHIC and LHC - the research programs 
concentrate on the study of the properties of deconfined QCD matter at very high 
temperatures and almost zero net baryon densities, whereas at moderate beam energies 
(SPS, FAIR and NICA) experiments focus on the search for structures in the QCD phase 
diagram such as the critical endpoint, the predicted first order phase transition between 
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Figure 1. The phase diagram with the critical end point at « 400 MeV,T « 
150 MeV as predicted by Lattice QCD. In addition, the time evolution in the 
T — /i-plane of a central cell in UrQMD calculations [MJ[3S] is depicted for different 
bombarding energies. At RHIC (see insert at the /z_b scale) large jiB can be accessible 
in the fragmentation region only. The figure is taken from Ref. [55] . 



hadronic and partonic matter, and the chiral phase transition. The critical endpoint and 
the first order phase transition are expected to occur at finite baryon chemical potential 
and moderate temperatures. 

Particle yields or ratios measured at different beam energies and analyzed with the 
statistical model provide sets of thermal parameters, temperature T and baryo-chemical 
potential /i^, which establish a "line of chemical freeze-out" [T6HT8]. The results of the 
fits to experimental data are shown in a phase diagram of hadronic and quark-gluon 
matter in Fig. [1] with full circles. 

Lattice QCD results [20l[2T] show a crossing, but no first order phase transition to 
the QGP for vanishing or small chemical potentials /is, i.e. at the conditions accessible 
at central rapidities at RHIC full energies. A first order phase transition is assumed 
at high baryochemical potentials or baryon densities. This corresponds to SPS energies 
and in the fragmentation region of RHIC |22||23]. The critical baryochemical potential is 
predicted [20l|2l] to be /i^ = 400 ± 50 MeV and the critical temperature = 140 ^ 160 
MeV. 

The thermodynamic parameters T and /i^ - as extracted from the microscopic 
UrQMD approach in the central overlap region of Pb+Pb or Au+Au collisions [26] - 
are shown in Fig. [T|, where the full dots with error bars denote the chemical freeze- 
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out parameters determined from fits to the experimental yields [TS]. The triangular 
and quadratic symbols (time-ordered in vertical sequence) stand for temperatures T 
and chemical potentials extracted from UrQMD transport calculations in central 
Au+Au (Pb+Pb) collisions at RHIC (21.3 ATeV), SPS (160, 40 AGeV) and AGS (11 
AGeV) p^l25] as a function of the reaction time (separated by 1 fm/c steps from top to 
bottom). The open symbols denote non-equilibrium configurations and correspond to 
T parameters extracted from the transverse momentum distributions, whereas the full 
symbols denote configurations in approximate pressure equilibrium in longitudinal and 
transverse direction. The T and /i^ values found in the transport model suggest that 
local thermal and chemical equilibrium may exist in the interesting regions of the phase 
diagram: the line of the first order phase transition, the critical point, and the region 
of a crossover. Note, however, that most transport models in their present formulation 
(such as UrQMD) have no explicit QGP phase. 

1.2. Fluctuations and Correlations in High Energy Collisions 

Fluctuations and correlations are important characteristics of any physical system. They 
provide essential information about the effective degrees of freedom. In general, one can 
distinguish between several classes of fluctuations. There are "dynamical" fluctuations 
and correlations reflecting the underlying dynamics of the system. There are also 
"trivial" fluctuations induced by the measurement process itself, such as finite number 
statistics, etc. They should be clarified and subtracted in order to access the dynamical 
fluctuations which tell as about the properties of the system as well as sucseptibilities. 

The event-by-event fluctuations in high energy A+A collisions are expected to be 
closely related to the transitions between different phases of QCD matter. By measuring 
the fluctuations one should observe anomalies from the onset of deconfinement and 
dynamical instabilities when the expanding system goes through the 1-st order transition 
line between the QGP and the hadron gas. It is predicted [27] that the onset of 
deconfinement should lead to a non-monotonous behavior in the energy dependence 
of multiplicity fluctuations, the so-called "shark fin". Furthermore, the QCD critical 
point may be signaled by a characteristic pattern in enhanced fluctuations |28j . 

Indeed, the Lattice QCD calculations in Ref. [29] indicate a sharp increase of the 
susceptibilities at the critical point for non-zero chemical potentials. This is illustrated 
in Fig. [2]which shows the dimensionless quark number susceptibility Xq/T"^ ^is a function 
of T/Tq for various Hq/T where the quark number (g) susceptibility is defined for the 
system in thermal equilibrium (for a grand-canonical ensemble) as follows [29] : 

Xq_ ^ ( d d \ n^ + UA , s 

T2 \^{^lJT)^ ^{^^JT)) T3 ^> 

with riq denoting the quark number density. The peak, which develops in Xq with 
increasing /x^, is a sign that fluctuations in the baryon density are growing when 
approaching the critical endpoint in the {fi,T). Thus A-l-A collisions in the SPS energy 
region are expected to be a suitable tool for a search of signatures for the critical point. 
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Figure 2. x?/^^ as a function of T/Tq for various /Xg/T according to Ref. [29] . 

An ambitious experimental program for the search of the QCD critical point has 
been started by the NA61 Collaboration at the CERN SPS [51130]. The program includes 
a variation in the atomic number A of the colliding nuclei as well as an energy scan. 
This allows to explore the phase diagram in the T — fiB plane near the critical point. 
One expects to 'locate' the position of the critical point by studying its 'fluctuation 
signals'. High statistics multiplicity fluctuation data will be taken for p+p, C+C, S+S, 
In+In, and Pb+Pb collisions at bombarding energies of Eiab=^0, 20, 30, 40, 80, and 
158 AGeV. 

The QCD critical point is expected to be experimentally seen as a non-monotonic 
dependence of the multiplicity fluctuations, i.e. a specific combination of atomic number 
A and bombarding energy Eiab could move the chemical freeze-out of the system close 
to the critical point and show a 'spike' in the multiplicity fluctuations. Since the 
available microscopic transport approaches Hadron-String-Dynamics (HSD) and Ultra- 
relativistic-Quantum-Molecular-Dynamics (UrQMD), which operate from lower SIS to 
top RHIC energies, do not include explicitly a phase transition from a hadronic to a 
partonic phase, we can not make a clear suggestion for the location of the critical point 
- this is beyond the scope of hadron-string models. However, this study should be 
helpful in the interpretation of the upcoming experimental data since it will allow to 
subtract simple dynamical and geometrical effects from the expected QGP signal. The 
deviations of the future experimental data from the HSD and UrQMD predictions may 
be considered as an indication for the critical point signals. 

We will deal in particular with particle number event-by-event fluctuations. They 
are quantifled by the ratio of the variance of the multiplicity distribution to its 
mean value, the scaled variance. Let's introduce the corresponding notations: The 
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deviation AA^^ from the average number (Na) of the particle species A is defined by 
Na = (Na) + AA^^, while the covariance for species A and B is 

A{Na,Nb) = {ANaANb) = {NaNb) - {Na){Nb) , (2) 

the variance, 

Var{NA) = A{Na,Na) = {{ANa)') = {ND - (Na)' , (3) 
the scaled variance, 

^ VarjNA) ^ {{ANAf) ^ {Nl) - {Na)' 
^ " {Na) {Na) {Na) ' ^ ^ 

and the correlation coefficient, 

{ANa ANb) 



Pab 



.2\ //A Ar x2\ll/2 • 



{{ANaY) {{ANbY) 



(5) 



Note that an absence of correlations, pab = 0, implies P{Na, Nb) = Pa{Na) x Pb{Nb), 
and w = 1 for the Poisson multiplicity distribution, P{N) = N exp (^— /N\ . 



1.3. Outline of the Review 

The Sec. [2] provides the basic concept of transport approaches. The centrality (Sec. [3]), 
energy (Sec. S]) and system size dependence (Sec. E]) of event-by-event fluctuations in 
the particle multiplicity are studied in A+A collisions at the SPS CERN energies. The 
results on particle number fluctuations and long-range forward-backward multiplicity 
correlations at RHIC BNL energies are discussed in Sec. El Electric charge and baryonic 
number fluctuations are considered in Sec. [7] while fluctuations of the particle number 
ratios are presented in Sec. [HI A summary in Sec. [9] concludes the review. 

The transport model calculations are performed within the HSD and UrQMD 
approaches. The results of both models qualitatively agree to each other. In some 
cases - to illustrate the systematic uncertainties of the transport model approaches for 
a speciflc observable - we will present the results of both models. Otherwise the results 
of only HSD simulations are presented for brevity. 



2. Concept of Transport Approaches 

Two independent transport models that employ hadronic and string degrees of freedom, 
i.e. HSD [311132] and UrQMD [33l|3l], have been used for studying the fluctuations in 
p+p and A+A collisions. Both approaches take into account the formation and multiple 
re-scattering of hadrons in A-l-A reactions and thus dynamically describe the generation 
of pressure in the hadronic expansion phase. Though quite different in the numerical 
realization, both models are based on the same concepts and dynamical degrees of 
freedom: strings, quarks, diquarks (g, g, gg, gg) as well as hadronic degrees of freedom 
in the conflned phase. We stress again that both approaches do not include an explicit 
phase transition to a quark-gluon plasma. 
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2.1. The covariant off-shell transport approach 

The novel HSD model (extended for off-shell dynamics [35]) is based on the solution 
of generalized covariant off-shell transport equations for hadronic degrees of freedom 
(for the details of the covariant off-shell transport theory we refer the reader to the 
review [36] and references therein). As demonstrated in [35] the off-shell dynamics is 
particularly important for the transport description of the propagation and interaction 
of resonances with broad spectral functions which change their properties in the hot and 
dense medium (e.g. for the vector mesons). 

The coupled set of transport equations for the phase-space distributions Nh{x,p) 
(x = it,r), p = {e,p)) of fermion h [3ll[32l|37] with a spectral function Ah{x,p) (using 
the Green function decomposition iG^{x^p) = Nh{x,p)Ah{x,p)) is formally written as 

{(n^ - u,d;u;: - mid^^u^) + {u,d;u;: + mid;u^) a^} n,{x,p) a,{x,p) 

- {zS<, Re G^} = Yl Tr^Tr.Tr, [T^T],2^,,6\U + Ha - - U,) (6) 

h2h3h4 

X Ah{x, p)Ah2 {x, P2)Ah^ {x, pzjAh^ (x, Pa) 

X {Nhz{x,p^)Nh4{x,Pi)fh{x,p)fhAx,P2) - Nh{x,p)Nh^{x,p2)fh-A^,P^)fhA^,Pi)] ■ 

Here = {dt, Vr), df^ = {d^, Vp), (/i = 0, 1, 2, 3) and fh{x,p) = 1-Nhix,p) for fermions. 

The 'backflow' term in ([6]) - i.e. the Poisson bracket with the self-energy and 
retarded Green function - is given by 

Mhix,p) 



p)2 + r4a;,p)74. 



d;(N,(x,p)r,Xx,p)) (7) 



-^K M,.(.Jtr:L)V4 )^'^'-'--''>^'-'--^>'- 

It stands for the off-shell evolution which vanishes in the on-shell limit, when the spectral 
function Ah{x,p) does not change its shape during the propagation through the medium, 
i.e. Vr.r(x,j9)=0 and VpT{x,p)=0 with T denoting the local width of the spectral 
function. 

In ([6]) the trace over particles 2,3,4 reads explicitly for fermions 
^2,r2 (27r)4 J 2Jpl + M| 



2 

where (T2, T2 stand for the spin and isospin of particle 2. In case of bosons one has 

2^^2 = E At / d'p2^, (9) 
(27r)4 y ^2 2 ' ^ ^ 

since here the spectral function Ab is normalized as 

-.2 



J^As{x,p) = l (10) 
whereas for fermions one obtains the normalization 

^A,{x,p) = l. (11) 
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For the bosons B the mass function is Eq. ([7]) 

Mb(x,p) =j9o-p ^-m^-ReS(x,p) (12) 

where the (local and nonlocal) self-energies are included in ReS(a;,p). 
The Lorentz covariant mass-function for the fermions is 

M,(x,p) = Ul-U'-mf , (13) 

with the effective mass and four-momentum given by 

mlix,p) = mh + U^{x,p) (14) 

U^{x,p) = p^-U}:{x,p) . 

In (HM U^{x,p) and Ul^{x,p) denote the real part of the scalar and vector hadron 
self-energies, respectively, and nih stands for the bare (vacuum) mass of baryon h. 

In ([6]) [T+T] 12^34 is the 'transition rate' for the process 1 + 2 — )■ 3+4 (which is taken 
to be on-shell in the default HSD approach). In the c.m.s. of the colliding particles the 
transition rate is given by 

da 



(15) 

1+2^3+4 

where da/dVt is the differential cross section of the reaction and Vi2 the relative velocity 
of particles 1 and 2. Note, that a generalized collision term for n ■H- m reactions is given 
in Ref. [38] and optionally included in HSD, which is, however, not used in the present 
study. 

The hadron quasi-particle properties in are defined via the mass-functions (fT2!) 
or (fT3l) with ( fT4l) while the phase-space factors 

Mx,p) = l±Nh{x,p) (16) 

are responsible for fermion Pauli-blocking or Bose enhancement, respectively, depending 
on the type of hadron in the final/initial channel. The transport approach ([6]) is 
fully specified by U^{x,p) and Ui^{x,p), which determine the mean-field propagation 
of the hadrons, and by the transition rates T'^T in the collision term, that describe 
the scattering and hadron production/absorption rates. The coupled set of hadronic 
transport equations may be employed for a large variety of hadronic systems at 
relativistic (and non-relativistic) energies. 



2.2. On-shell limit 

In our present study we will discard the scalar and vector potentials and Uj^ in ([6]) 
and (fT^ which are of minor importance at SPS and RHIC energies due to the dominance 
of string excitations and decay. Furthermore, we will employ the on-shell limit for the 
spectral functions 

Ah{x,p) r^6{p'^ - ml), 

where ruh denotes the hadron on-shell mass. These (legitimate) approximations - 
employed in the standard transport approximations - will allow for a large reduction 
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in the necessary CPU time and a proper statistics for the fluctuation analysis. We 
mention that each individual cascade simulation will correspond to a micro-canonical 
non-equilibrium simulation of the dynamics with exact conservation of energy and 
momentum, baryon number, charge number, strangeness etc. 

Also note that the transport equations - in the limits discussed above - yield the 
well-known statistical distributions 



with ?7 = 1 for mesons and rj = —1 for fermions in the limit t — )• oo, if the 
system is conflned to a flnite (and flxed) volume. In (fTTll fih and T denote the 
conventional Lagrange parameters for the chemical potential of hadron h and equilibrium 
temperature. Accordingly, the transport approaches allow to explore non-equilibrium 
dynamical effects for fluctuations observables while they smoothly approach the grand 
canonical ensemble (GCE) in the limit t — > oo for stationary systems when averaging 
over many events with a different amount of excitation energy and different number of 
produced particles. 

2.3. Basic concept of HSD and UrQMD 

We recall that in the HSD (v. 2.5) [3ll[32l[39] approach nucleons, A's, N*(1440), 
N*(1535), A, S and S* hyperons, H's, S*'s and fi's as well as their antiparticles are 
included on the baryonic side whereas the 0~ and 1~ octet states are incorporated in 
the mesonic sector. Inelastic baryon-baryon (and meson-baryon) collisions with energies 
above ^/s^|^ ~ 2.6GeV (and ^/sth — 2.3GeV) are described by the FRITIOF string model 
|40] whereas low energy hadron-hadron collisions are modelled in line with experimental 
cross sections. Low energy cross sections such as threshold meson production in proton- 
neutron (pn) collisions - which are scarcely available from experiments - are flxed 
by proton-proton (pp) cross sections and isospin factors emerging from pion-exchange 
diagrams. Since we address ultra-relativistic collisions at SPS and RHIC energies 
such 'low energy uncertainties' are of minor relevance here. As pre-hadronic degrees 
of freedom HSD includes 'effective' quarks (antiquarks) and diquarks (antidiquarks) 
which interact with cross sections in accordance with the constituent quark model (cf. 



The UrQMD (v. 1.3) transport approach [331IM] includes all baryonic resonances up 
to masses of 2 GeV as well as mesonic resonances up to 1.9 GeV as tabulated by the 
Particle Data Group [12]. For hadronic continuum excitations a string model is used 
with hadron formation times in the order of 1-2 fm/c depending on the momentum and 
energy of the created hadron. The transport approach - by construction - yields a good 
reproduction of nucleon-nucleon, meson-nucleon and meson-meson cross section data in 
a wide kinematical range [33|IM]. 

The HSD and UrQMD transport models have been used for the description of pp, 
pA and AA collisions from SIS to RHIC energies and lead to a fair reproduction of 'bulk' 
properties such as hadron abundances, rapidity distributions etc. We mention that some 




(17) 



Refs. |4T]). 
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rare probes (e.g. particles with high transverse momentum) are not well described by 
transport models (cf. Refs. |26lll3llll]), however, this issue is of minor impact on the 
fluctuation analysis to be presented below. 

3. Multiplicity Fluctuations in Nucleus-Nucleus Collisions at the SPS 

This Section contains the results for the particle number fluctuations in Pb+Pb collisions 
at 158 AGeV within HSD and UrQMD |15]. The results of these transport models are 
compared with the model of independent sources and with the experimental data of the 
NA49 Collaboration. 

3.1. Fluctuations in the Number of Participants 

Fig. |3] demonstrates a typical non-central A+A collision. In each collision only a 
fraction of all 2 A nucleons interact. These are called participant nucleons and are 
denoted as N^"-' and N^p"^^ for the projectile and target nuclei, respectively. The 
corresponding nucleons, which do not interact, are called the projectile and target 
spectators, N^"-' = A — N'p°-' and A^^"^^ = A — N^p^^ . In the transport models the 
nucleon spectators are deflned according to the criteria, \y — ybeam(target)\ < 0.32, which 
restricts their rapidities. 




Figure 3. Non-central heavy ion collision: Spectators in the projectile and target are 
denoted as Np^"-' and N^p^^ . 'Participant region' indicates the region where the particle 
production occurs due to the interaction of participants from target and projectile. 
Note that in a fixed target experiment it is possible to fix the number of projectile 
participants by detecting projectile spectators via a Zero Degree Calorimeter (ZDC). 

The fluctuations in high energy A+A collisions are dominated by a geometrical 
variation of the impact parameter. However, even for the flxed impact parameter the 
number of participants, Np = Np "-' + N'^^^ , fluctuates from event to event. This is 
due to the fluctuations of the initial states of the colliding nuclei and the probabilistic 
character of the interaction process. The fluctuations of Np form usually a large and 
uninteresting background. In order to minimize its contribution the NA49 Collaboration 
has selected samples of collisions with a flxed numbers of projectile participants. This 
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selection is possible due to a measurement of N'^'"-' in each individual A+A collision by 
a calorimeter which covers the projectile fragmentation domain. However, even in the 
samples with N'p"-' = const the number of target participants fluctuates considerably. 
Hence, an asymmetry between projectile and target participants is introduced, i.e. Np'"-' 
is constant by constraint, whereas N^p^^ fluctuates. 







Pb+Pb, 158 A GeV p 



HSD 
UrQMD 




50 



100 



150 



200 



N 



proj 



Figure 4. The HSD and UrQMD results for the scaled variance uj*p^^ as the function 



From an output of the HSD minimum bias simulations of Pb+Pb collisions at 
158 AGeV the samples of events with fixed values of Np'"-' have been formed. For the 
average number of target participants in these samples one finds, {N^p^^) = N^"-' . The 
deviations are only seen at very small {Np'"-' ^ 1) and very large {Np'"-' ^ A) numbers 
of projectile participants. In each sample with fixed value of Np'"-' the N'^^^ number 
fluctuates around its mean value. The scaled variance of these fluctuations is denoted 
as uip^^. Fig. m shows both the HSD and UrQMD results for the scaled variance Up^^ 
as the function of Np""-^ . The fluctuations of Np^^ are quite strong; the largest value of 



targ 



3 ^ 3.5 occurs at Nl^"^ = 20 ^ 30. 



3.2. HSD and UrQMD Results for Pb+Pb at 158 AGeV 



The NA49 Collaboration has minimized the event by event fluctuations of the number of 
nucleon participants in measuring the multiplicity fluctuations by selecting the samples 
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of collisions with a fixed number of projectile spectators, A'^^'^"-' = const, and thus a 
fixed number of projectile participants, N'p"-' = A — N'g'"-' . From an output of the HSD 
and UrQMD minimum bias simulations the samples of Pb+Pb events with fixed values 
of Np^"-' have been formed. Fig. [5] presents the HSD and UrQMD results for the scaled 
variances of negatively, positively, and all charged particles, Ui (with i = +,—,ch), in 
Pb+Pb collisions at 158 AGeV calculated at fixed Np'"-' . 

The final particles in the HSD and UrQMD simulations are accepted at rapidities 
1.1 < y < 2.6 {y corresponds to the particle rapidity in the Pb+Pb c.m.s. frame) in 
accordance to the NA49 transverse momentum filter [46]. This is done to compare the 
HSD and UrQMD results with the NA49 data. The HSD and UrQMD simulations both 
show fiat Ui values, ~ a;+ ~ 1.2, Uch ~ 1-5, and exhibit almost no dependence on 
N^"-' . The NA49 data, in contrast, exhibit an enhancement in Ui for N^"-' ^ 50. The 
data show maximum values, ~ w+ ~ 2 and Uch ~ 3, and a rather strong dependence 
on NP;'"^. 

Fig. [5] also shows results of the HSD and UrQMD simulations for the full Air 
acceptance for final particles, and shows the NA49-like acceptance in the mirror rapidity 
interval, —2.6 < y < —1.1 of the target hemisphere. HSD and UrQMD both result 
in large values of Ui, i.e. large fluctuations in the backward hemisphere: in the 
backward rapidity interval —2.6 < y < —1.1 (target hemisphere) the fluctuations 
are much larger than those calculated in the forward rapidity interval 1.1 < y < 2.6 
(projectile hemisphere, where the NA49 measurements have been done). Even much 
larger fluctuations seen in Fig. [5] follow from the HSD and UrQMD simulations for the 
full acceptance of final particles. 

At fixed values of the numbers of participants Np^"-' and Np^^ one can introduce 
the probability Wi{Ni; Np^^ , Np"-') for producing Ni {i = —,+,ch) final hadrons. At 
fixed N^°-^ the averaging procedure is defined as, 

(■■■)= E E ■ ■ ■ WpiN';''; N^') W,{N,; N'^^^ , N^') , (18) 

where Wp{N^p^^; Np'"-') is the probability for a given value of N'^^^ in a sample of events 
with fixed number of the projectile participants Np'"-' . Under certain simplifications, 
the variance can be presented as, 

Var{N,) = {N^) - {N,f = u* {N,) + up m (Ni) , (19) 

where the scaled variance u* corresponds to the fluctuations of at fixed values of 
N^p'"' and iVp'^^ m = {Ni)/{Np) , and Np = N'^i'^ + N"";"^ is the total number of 
participants. In obtaining of Eq. ( IT9|) two assumptions have been made. First, it is 
assumed that u* does not depend on Np. The second assumption is that the average 
multiplicities {Ni) are proportional to the number of participating nucleons. Finally, 
the scaled variances, cjj, can be presented as: 

VariN.^ , , , 

OJi = ^ = + upUi . (20) 
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Figure 5. The results of the HSD () and UrQMD (right) simulations are shown for 
a;_, and LOch in Pb+Pb collisions at 158 AGeV as functions of N^°-' . The black 
points are the NA49 data. The different lines correspond to the model simulations 
with the original NA49 acceptance, 1.1 < y < 2.6, in the projectile hemisphere (lower 
lines), the NA49-like acceptance in the mirror rapidity interval, —2.6 < y < —1-1, in 
the target hemisphere (middle lines), and in full An acceptance (upper lines). 
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The average values are {N^p''^) = N'p"-' , thus, {N^p"^^) = {Np)/2. Therefore, the scaled 
variance up for the total number of participants in Eq. fl20|) equals to up = (Jp"^^ /2 
as only a half of the total number of participants does fluctuate. The value of u/p^^ 
depends on Np"-' , as shown by the HSD and UrQMD results in Fig. HI The average 
particle number = +, — , ch) per participant calculated within the HSD and UrQMD 
models for full 47r-acceptance show a good agreement with each other as well as with 
the extrapolated to 47r NA49 data. 

Eq. ( I20i) corresponds to the so called model of independent sources (MIS), see, e.g., 
Ref [U]. One assumes that final particles are produced by the sources. The numbers 
of the sources are taken to be proportional to the number of projectile and target 
participant nucleons, and the sources are assumed to be independent of each other. 
The physical meaning of the sources depends on the model under consideration (e.g., 
wounded nucleons [48j, strings and resonances [STtiM] . or the fluid cells at chemical 
freeze-out in the hydrodynamical models). The Eq. (120|) presents the final multiplicity 
fluctuations as a sum of two terms: the fluctuations from one source, cj*, and the 
contribution due to the fluctuations of the number of sources, upiii. 

In peripheral A+A collisions there are only few nucleon-nucleon (N+N) collisions, 
and rescatterings are rare. The picture of independent N+N collisions looks thus 
reasonable. In this hadron production source can be associated with a N+N 

collision. The multiplicities and scaled variances of N+N is given in terms of proton- 
proton (p+p), proton-neutron (p+n), and neutron- neutron (n+n) collisions: 

(AT/^^) = a,, {Nf) + (ATf") + a„„ {ND , (21) 

^ = [o^vv {Nf) + ujT {ND + «„„ ur {ND] , (22) 

where a^p = Z^/A^ = 0.155, a^n = 2Z{A~Z)/A^ = 0.478, a„„ = {A-Zf/A^ = 0.367 
are the probabilities of p+p, p+n, and n+n collisions in Pb+Pb reactions (A=208, 
Z=82). The average multiplicities and scaled variances for elementary collisions 
calculated within the HSD simulations at 158 GeV are equal to: 

(AT^n = 6.2 , {ND = 5.8 , {ND = 5.4 , (23) 

coZ = 2.1 , cor, = 2.4 , = 2.9 . (24) 

For negatively and positively charged hadrons, the average multiplicities and scaled 
variances in elementary reactions can be presented in terms of corresponding quantities 
for all charged particles: 

(iV±) = I {N.h) ± 7) , 00^ = I ^cH , (25) 

2 2 {Nch} ± 7 



with 7 = 2, 1, for pp, pn and nn reactions, respectively. Thus, using Eqs. fl2m25p one 
finds the HSD results for u* per N+N collision at 158 GeV: 

uj*f^ = 2.5 , = 1.5 , ul = 1.1 . (26) 

The above arguments are not applicable for central A+A collisions, where a large 
degree of thermalization is expected. In the limit of Np"-' = A one can take the values 
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of u* from the Pb+Pb data or model simulations. In this limit, ujp = Up^^ /2 ^ (see 
Fig. H]), and thus Ui ~ u*. It has been found that fl26l) gives a reasonable description 
of Ui in the HSD simulations for central Pb+Pb collisions at 158 AGeV. Therefore, one 
can extrapolate (!20|) and ( 126|) to all values of Np"-' . This leads to a good agreement for 
Ui (l20l) in the full Att acceptance with the transport model simulations [15]. 

The fluctuations of the total hadron multiplicities generated by the transport model 
dynamics are large, i.e. the Ui for 47r-acceptance are essentially larger than 1. The main 
contributions to Ui come from the second terms in ( l20i) , which are due to the fluctuations 
of N'^^^. These fluctuations of the target nucleon participants presented in Fig.lHexplain 
both, the large values of Ui and their strong dependence on Np'"-' . 



4 
3 

3 2 



1 


50 100 150 200 

lyproj ^ -^targ 




Figure 6. The circles, triangles, and boxes are the results of the HSD simulations for 
uji in full 47r acceptance with N*p^^ = Np'"-' . This condition yields Lu*p^^ = 0, and 
Eq. (|20|) is reduced to uJi — uj* . The dashed lines correspond to w* taken from Eq. (|26|) . 

Fig. |6] supports the previous findings. The HSD events with fixed target participant 
number, N^^^ = N^°\ exhibit much smaller multiplicity fluctuations. This is due to the 
fact that terms proportional to u/p^^ in (120|) do not contribute, Ui become approximately 
independent of the number of participants and equal to uj* (l26ll . 

3.3. Transparency and Mixing in Nucleus-Nucleus Collisions 

Different models of hadron production in relativistic A+A collisions can be divided into 
three limiting groups: transparency, mixing, and reflection models (see Ref. [49]). The 
flrst group assumes that the flnal longitudinal flows of the hadron production sources 
related to projectile and target participants follow in the directions of the projectile and 
target, respectively. One calls this group of models as transparency (T-)models. If the 
projectile and target flows of hadron production sources are mixed, these models are 
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y 
target projectile 

Figure 7. The sketch of the rapidity distributions of the baryon number or the particle 
production sources (horizontal rectangles) in nucleus-nucleus collisions resulting from 
the transparency, mixing and reflection models. The spectator nucleons are indicated 
by the vertical rectangles. In the collisions with a fixed number of projectile spectators 
only matter related to the target shows significant fluctuations (vertical arrows). See 
Ref. [49] for more details. 

called as mixing (M-)models. Finally, one may assume that the initial flows are reflected 
in the collision process. The projectile related matter then flows in the direction of the 
target and the target related matter flows in the direction of the projectile. This class 
of models corresponds to the reflection (R-)models. The rapidity distributions resulting 
from the T-, M-, and R-models are sketched in Fig. [7] taken from Ref. ^9] . 

An asymmetry between the projectile and target participants introduced by the 
experimental selection procedure in a flx target experiment can be used to distinguish 
between projectile related and target related flnal state flows of hadron production 
sources. The multiplicity fluctuations measured in the target momentum hemisphere 
clearly are larger than those measured in the projectile hemisphere in T-models. The 
opposite relation is predicted for R-models, whereas for M-models the fluctuations 
in the projectile and target hemispheres are expected to be the same. Note that 
there are models which assume the mixing of hadron production sources, however, the 
transparency of baryon flows, e.g. three- fluid hydrodynamical model [501151]. 

Now the fluctuations of the particle multiplicities in the projectile {y > 0) and 
target {y < 0) hemispheres are considered. As one can see from Fig. HI in samples with 
Np'"-' = const the number of target participants, N^p^^ , fluctuates considerably. Of 
course, this event selection procedure introduces an asymmetry between projectile and 
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Figure 8. The scaled variances oji for the projectile (boxes) and target (circles) 
hemispheres in the HSD {left) and UrQMD {right) simulations. 

target participants: N'p°-' is constant, whereas N^p'^^ fluctuates. Under this restriction 
the transport models give very different results for the particle number fluctuations in 
the projectile and target hemispheres. As clearly seen from Fig. [8] the particle number 
fluctuations in the target hemispheres are much stronger than those in the projectile 
hemispheres. There is also a strong iVp^""'- dependence of Ui in the target hemisphere, 
which is almost absent for the ui in the projectile hemisphere. Thus, the fluctuations 
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Figure 9. Scaled variances a;_,cj+, and uJch from ttie HSD (left) and UrQMD (riglit) 
in the rapidity intervals Ay = 1 around y = 0, ±0.5, ±1.5, ±2.5 for different number of 
projectile participants N^"-^*. 



of N'^^^ have a small influence on the final multiplicity fluctuations in the projectile 
hemisphere, but they contribute very strongly to those in the target hemisphere. 

A selection of collisions with a fixed number of Np'"-' and fluctuating number of 
N'^^^ means that the projectile and target initial flows are marked in fluctuations [19] in 
the number of colliding nucleons. The projectile and target related matters in the final 
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state of collisions can be then distinguished by an analysis of fluctuations of extensive 
quantities. The analysis of the fluctuations applied to collisions of identical nuclei gives 
a unique possibility to investigate the flows of particle production sources. 

Fig. [9] shows the scaled variances u_,(jJ+ and cOch in the HSD and UrQMD 
simulations as functions of rapidity for different N^"^"^ values. It is clearly seen that 
the bias on a fixed number of projectile participants reduces strongly the particle 
fluctuations in the forward hemisphere, in particular within the NA49 acceptance 
(1.1 < y < 2.6). The fluctuations of the target participant numbers influence strongly 
the hadron production sources in the target hemispheres. They also contribute to the 
projectile hemisphere, but this contribution is only important in the rapidity interval 

< y < 1, i.e. close to midrapidity. It turns out that this "correlation length" in 
rapidity. Ay ^ 1, a.s seen in Fig. [H is not large enough to reproduce the data. The 
large values of Ui and their strong A^'p'^"-'- dependence in the NA49 data (cf. Fig. [5]) 
in the projectile rapidity interval, 1.1 < y < 2.6, thus demonstrate a significantly 
larger amount of mixing in peripheral reactions than generated in simple hadron/string 
transport approaches. 

4. Multiplicity Fluctuations in Central Nucleus-Nucleus Collisions 

This Section presents the HSD results on the excitation function of the multiplicity 
fluctuations in central A+A collisions. They will be compared with corresponding results 
in p+p collisions. Note, that some observables in p+p and A+A collisions are rather 
close to each other. For example, the charged hadron multiplicity per participating 
nucleon, Uch = (Nch) / (^p)? at SPS energies of 20 -r- 158 AGeV are not much different in 
central Pb+Pb and inelastic p+p collisions [l3l[52], Rch = {^chjAAl i^ch)pp = 1 ^ 1.5 0. 
This explains a vitality of the wounded nucleon model (WNM) |48j which treats 
the final state in A+A collision as the result of independent N+N collisions (see, 
e.g., Ref. [5l] which discusses the recent data [55] on d+Au collisions at the RHIC 
energies of ^/sNN = 200 GeV). However, the basic concept of the WNM is in a 
severe conflict with many other data, e.g., with multi-strange baryon production, 
Rn = (^q)aa / (^^^)pp - 12.5 [56] in Pb+Pb at 158 AGeV. The Rn enhancement 
is expected to be even stronger at smaller collision energies. The search for quark- 
gluon plasma signatures in A+A collisions is usually based on the expectation of a very 
different behavior of special physical observables in A+A and p+p collisions. Famous 
examples of QGP signatures are the 'strangeness enhancement', 'J/ip suppression', and 
'jet quenching'. In all these cases one compares a suitably normalized physical quantity 
in A+A and in p+p reactions at the same collision energy per nucleon. 

In general, one can define two groups of hadron observables. The first group includes 
observables which are rather similar in A+A and p+p collisions, thus, they can be 
reasonably described within the WNM. The second group consists of A+A observables 
which are very different from those in p+p collisions. The question arises: are the 

1 Note that Rch < 1 at low collision energies. The change in behavior of Rch is discussed in Ref. [53] . 



4- Multiplicity Fluctuations in Central Nucleus- Nucleus Collisions 



22 



multiplicity fluctuations in A+A collisions close to those in p+p reactions, or are they 
very different? The aim of present Section is to study this question [57] . 



4.1. Inelastic N+N and Central A+A Collisions 

The compilation of p+p data for (Nch) and Uch are taken from Ref. [17] and presented 
in Fig. [101 The energy dependence can be parameterized by the functions [17] : 

0.31 /,Ar \ 



-4.2 . 4.9 (^) 



0.35 



mn) - If 



(27) 



At high collision energies the KNO scaling [58] holds which implies that the multiplicity 
distribution P{Nch) behaves as 

(see also Ref. [591160]). For (Nch) > 1 it follows, 

In particular, Uch (Nch) [61] as also seen from the parametrization ( 1271) . 
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Figure 10. The multiplicity (left) and scaled variance (right) of all charged hadrons 
in p+p inelastic collisions as functions of collision energy. The dashed lines correspond 
to the parametrization ((27l) from Ref. [471. The solid lines are the HSD results. 

The HSD model description of the p+p data (for p+p reaction this is almost 
equivalent to the Lund-String model [62]) is shown in Fig. [TOl by the solid lines. It 
gives a good reproduction of the p+p data for {Nch), but slightly underestimates Uch 
at high collision energies. For negatively and positively charged hadrons the average 
multiplicities and scaled variances in p+p collisions can be presented in terms of the 
corresponding quantities for all charged particles, see Eq. fl25l) . 

We compare now the central collisions of heavy nuclei and N+N inelastic collisions 
( I2TII22I) within the HSD model. A small difference between p+p and N+N collisions is 
only present at the SPS energies and gradually disappears at RHIC energies. 
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Figure 11. The multiplicities per participant, Ui (left), and scaled variances, u)i 
(right). The solid lines are the HSD results for N+N collisions according to (I^Tt . The 
full circles are the HSD results for central A+A collisions for zero impact parameter, 
b = 0. The full squares for n_ are the NA49 data [l3l[52] for ((tt") + {K-))/{Np) in 
the samples of 7% most central Pb+Pb collisions. The HSD results for oj^ after the 
subtraction of the contributions Aw^ = riiUjp. are shown by open triangles. The dotted 
lines are the MCE HG model results for Wi [63]. The HG parameters correspond to 
the chemical freeze-out conditions found from fitting the hadron yields. 

In Fig. [11] the HSD model results are shown for the multiplicities per participating 
nucleons, rij = {Ni)/{Np), and for the scaled variances, Ui, in central collisions (zero 
impact parameter, 6 = 0) of Pb+Pb at Eiab=lO, 20, 30, 40, 80, 158 AGeV and Au+Au at 
ysjv7v=62, 130, 200 GeV. From Fig.[TT]one concludes that the HSD results for the scaled 
variances in central A+A collisions are close to those in inelastic N+N collisions. The 
participant number fluctuations are found to be rather small for central A+A collisions 
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with zero impact parameter 6 = 0. For example, in Pb+Pb collisions with 6 = at 
158 AGeV the mean number of participants is (Np) = 392, and the scaled variance 
is cjp = 0.055 . The additional fluctuations, Aui, of ith hadrons due to participant 
number fluctuations can be estimated as, Aui = Ui ujp . The HSD results for uji after 
subtraction of the contributions Acjj are shown in Fig. [TT] by open triangles. These 
contributions to uji due to participant number fluctuations appear to be rather small 
at 6 = 0, and they do not explain the (positive) difference, a;j(AA) — a;j(NN) seen in 
Fig. [n] at = 200 GeV. 

On the other hand in the statistical model the scaled variances Ui = 1 for the 
ideal Boltzmann gas in the grand canonical ensemble (GCE). The deviations of Ui from 
unity in the hadron-resonance gas (HQ) model stem from Bose and Fermi statistics, 
resonance decays, and exactly enforced conservations laws within the canonical ensemble 
(CE) or micro-canonical ensemble (MCE) jSSHBB] . Note that the statistical model 
gives no predictions for the energy dependence of hadron multiplicities. All yields are 
proportional to the system volume V which is a free model parameter fitted to the 
multiplicity data at each collision energy. However, the statistical model does predict 
the scaled variances Ui as they become to be independent of the system volume for large 
systems. In Fig. [11] the scaled variances Ui calculated within the MCE HG model along 
the chemical freeze-out line (see Ref. [63] for details) are presented by the dotted lines: 
oji reach their asymptotic values at RHIC energies, a;±(MCE)= 0.3 and a;c/i(MCE)= 0.6. 
The corresponding results in the GCE and CE are the following: C(;±(GCE)= 1.2 and 
Wch(GCE)= 1.6, a;±(CE)= 0.8 and a;ch(CE)= 1.6. The HSD results for Ui in central 
A-l-A collisions are very different. They remain close to the corresponding values in 
p+p collisions and, thus, increase with collision energy as Ui oc Ui. One observes no 
indication for 'thermalization' of fluctuations in the HSD results. This is especially seen 
for RHIC energies: a;,(HSD)/wi(MCE)> 10 at ^JJ^ = 200 GeV. 

4-2. Comparison with the NA49 Data 

The fluctuations of the number of nucleon participants correspond to volume 
fluctuations, hence, they translate directly to the final multiplicity fluctuations. To avoid 
these 'trivial' fluctuations, one has to select a sample of very central, < 1%, collisions. 
Such a rigid centrality selection has been recently done for the NA49 data |66] by fixing 
the number of projectile participants close to its maximal value Np^"-' = A. 

The HG model was compared with the NA49 data [66] for the sample of 1% most 
central collisions at the SPS energies, 20 -r- 158 AGeV in Ref. [63]. It was found that 
the MCE results for u± are very close to the data, they are shown by the dashed lines 
in Fig. [121 The NA49 acceptance probabilities for positively and negatively charged 
hadrons are approximately equal, and their numerical values are: q = 0.038, 0.063, 
0.085, 0.131, 0.163, at the SPS energies of 20, 30, 40, 80, 158 AGeV, respectively. In 
the statistical model the scaled variances u'^'^ for the accepted particles are calculated 
from u± in the full space according to the acceptance scaling formulae (see Ref. ^3] for 
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details) : 

ul^" = 1 - q + q uj± . (28) 

Note that the energy dependence of u'^?'^ seen in Fig. [12] is strongly influenced by an 
increase with energy of the acceptance parameter q: only about 4% of the hadrons are 
detected at 20 AGeV and 16% at 158 AGeV. 
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Figure 12. Upper panel. The scaled variances u'^'' for central Pb+Pb collisions. The 
squares with error bars are the NA49 data for 1% most central collisions [66]. The 
dotted lines show the MCE HQ model results calculated from full Air scaled variances 
using ([28]) . The full circles present the HSD results in Pb+Pb collisions for 6 = with 
the NA49 experimental acceptance conditions, while the open circles are obtained from 
the 47r HSD scaled variances using ([28l) . Lower panel. The MCE HG (dotted line) and 
HSD (full circles, the same as in Fig. [TT]) for the Att scaled variances uj± at the SPS 
energies. 



The comparison of the HSD results for central Pb+Pb collisions (zero impact 
parameter, 6 = 0) with the NA49 data of 1% most central collisions, selected by 
the number of projectile spectators, is presented in Fig. [12] It demonstrates a good 
agreement of the HSD results with the NA49 data. There are also no essential differences 
between the MCE HG model and the HSD transport model results. Several comments 
are needed at this point. The HSD results within the NA49 acceptance demonstrate 
that the acceptance scaling formulae ([28|l is violated. The straightforward calculations 
(full circles in Fig. [T2]) lead to smaller values of u'^?'^ than those obtained with the 
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acceptance scaling formulae f l28p (open circles in Fig. [T^ . This difference may lead to 
a 10% effect in u'^'^ for the NA49 acceptance conditions. Thus, the MCE results for 
oj'j?'^ may also be about 10% smaller than those obtained from ( l28l) and shown in the 
upper panel of Fig. [121 The lower panel of Fig. [12] demonstrates that the MCE and 
HSD results for u± at the lowest SPS energy 20 AGeV are 'occasionally' rather close to 
each other. They both are also close to u± in p+p collisions. The HSD scaled variances 
Ui increase with collision energy. In contrast, the MCE cuj values remain approximately 
constant. The ratio of the HSD to MCE values of u± reaches about the factor of 2 at 
the highest SPS energy 158 AGeV. It becomes a factor of 10 at the top RHIC energy 
^/SNN = 200 GeV. However, the rigid centrahty selection is absent for the available 
RHIC fluctuation data. Due to this reason the participant number fluctuations give a 
dominant contribution to Wj. On the other hand, for the SPS data the small values of 
the acceptance, q = 0.04 -=-0.16, and 10% possible ambiguities coming from f l25]) almost 
mask the difference between the HSD and MCE results (Fig. [121 upper panel). 

Note that the scaled variance is a non-trivial function of the selected phase-space. 
This issue has been addressed in Ref. [67]. In order to study the dependence of scaled 
variance on rapidity, 12 different rapidity intervals have been constructed in such a way 
that the mean multiplicity in each interval is the same. If the scaled variance would 
follow the acceptance scaling formula ( 128]) . the scaled variance would be the same in 
each interval. However, as it has been shown in Ref. |67], the scaled variance is much 
higher near midrapidity than in forward and backward rapidities, which is illustrated in 
Fig. [T3] taken from Ref. [67J . 



~ 1.2| 1 1.2| 1 — 1.4| 

"3 • UrQMD, b=0 3 • UrQMD, b=0 "3 • UrQMD, b=0 



•••• 



Figure 13. Rapidity dependence of scaled variance in UrQMD simulation performed 
in full acceptance of positive (top), negative (middle) and all charged (bottom) hadrons 
in central Pb+Pb collisions at 158 AGeV. The rapidity bins are constructed in such 
a way that the mean multiplicity in each bin is the same. The figure is taken from 
Ref. |67|. 



The transverse momentum dependence of scaled variance is shown in Fig.[Tl]for the 
full longitudinal phase-space and for a midrapidity and a forward rapidity interval. The 
scaled variance decreases with increasing transverse momentum for the full acceptance 
and at forward rapidity. At midrapidity it stays approximately constant. The decrease 



5. Dependence on Energy and Atomic Number 



27 



"5" • 4lt 3 •41 "5" • 4ll 

1-1" A 0<y<0.5 1-1" A 0<y<0.5 A 0<;y<0.5 

T 1.25<y<1.75 T 1.25<y<1.75 ''■^^ T 1.25<y<1.75 



1 1.5 0.5 1 1.5 0.5 1 1.5 

(GeV/c) (GeV/c) p^ (GeV/c) 



Figure 14. (Color online) Transverse momentum dependence of multiplicity 
fluctuations of positively (top) , negatively (middle) and all charged hadrons (bottom) 
for all rapidities), < y < 0.5 and 1.25 < y < 1.75. The transverse momentum bins 
are constructed in such a way that the mean multiplicity in each bin is the same. The 
flgure is taken from Ref. [67]. 



of scaled variance is stronger for positively charged hadrons than for negatively charged 
ones because the protons, which have smaller relative fluctuations due to the large 
number of protons which enter the collision, have a larger mean transverse momentum. 

A similar effect of decreasing fluctuations for larger rapidities and transverse 
momenta is observed as a result of energy- and momentum conservation in a hadron gas 
model using the micro-canonical ensemble [M]. It costs more energy to create a particle 
with high momentum, therefore their number is expected to fluctuate less. 

5. Dependence on Energy and Atomic Number 

An ambitious experimental program for a search of the QCD critical point has been 
started by the NA61 Collaboration at the SPS |5|l30]. The program includes a variation 
in the atomic mass number A of the colliding nuclei as well as an energy scan. This 
allows to scan the phase diagram in the plane of temperature T and baryon chemical 
potential fiB near the critical point as argued in Ref. [5l[30] and shown in Fig. [15] which 
is taken from Ref. [69]. One expects to 'locate' the position of the critical point by 
studying its 'fluctuation signals'. High statistics multiplicity fluctuation data will be 
taken for p+p, C+C, S-l-S, In-|-In, and Pb-|-Pb collisions at bombarding energies of 
Eiab=lO, 20, 30, 40, 80, and 158 AGeV. 

The aim of the Section is to present the results of studies on the energy and system 
size dependence of event-by-event multiplicity fluctuations within the HSD and UrQMD 
microscopic transport approaches. Our study thus is in full correspondence to the 
experimental program of the NA61 Collaboration [5|[30]. 

The QCD critical point is expected to be experimentally seen as a non-monotonic 
dependence of the multiplicity fluctuations, i.e. a specific combination of atomic mass 
number A and bombarding energy Eiab could move the chemical freeze-out of the system 
close to the critical point and show a 'spike' in the multiplicity fluctuations. Since HSD 
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Figure 15. Left: The data sets on central A+A collisions planned to be registered 
by NA61 in a search for the critical point of strongly interacting matter and a study 
of the properties of the onset of deconfinement. Right: Hypothetical positions of the 
chemical freeze-out points of the reactions (In+In, S+S, C+C and p+p from bottom 
to top at 158A, 80A, AOA, 30A, 20A and lOA GeV from left to right) to be studied by 
NA61 in the (temperature)- (baryon-chemical potential) plane are shown by full dots. 
The open squares show the existing Pb-|-Pb NA49 data. 



and UrQMD do not include explicitly a phase transition from a hadronic to a partonic 
phase, a clear suggestion for the location of the critical point can not been made - 
it is beyond the scope of such hadron-string models. However, this study might be 
helpful in the interpretation of the upcoming experimental data since it will allow to 
subtract simple dynamical and geometrical effects from the expected QGP signal. The 
deviations of the future experimental data from the HSD and UrQMD predictions may 
be considered as an indication for the critical point signals. 

Theoretical estimates give about 10% increase of the multiplicity fluctuations due 
to the critical point [28 1I701I7T] . It is large enough to be observed experimentally within 
the statistics of NA61 [5l[30]. To achieve this goal, it is necessary to have a control 
on other possible sources of fluctuations. One of such sources is the fluctuation of the 
number of nucleon participants. It has been shown in Section [3] that these fluctuations 
give a dominant contribution to hadron multiplicity fluctuations in A-l-A collisions. 
It was demonstrated that one can suppress the participant number fluctuations by 
selecting most central A+A collisions. That's why the NA61 Collaboration plans to 
measure central collisions of light and intermediate ions instead of peripheral Pb-f-Pb 
collisions. It is important to stress, that the conditions for the centrality selection in the 
measurement of fluctuations are much more stringent than those for mean multiplicity 
measurements. 

5.1. Why Does One Need Central Collisions of Light Ions? 

To minimize the event-by-event fluctuations of the number of nucleon participants in 
measuring the multiplicity fluctuations the NA49 Collaboration has been trying to fix 
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Np'°-' in Pb+Pb collisions. Samples of collisions with a fixed number of projectile 
spectators, A^|^°-' = const (and thus a fixed number of projectile participants, Np'°-'), 
have been selected. A similar centrality selection is expected to be implemented in the 
NA61 experiment. 




Figure 16. The scaled variance uj*p^^ for tlic fluctuations of the number of target 
participants, N*p^^ . The HSD simulations of ajp*^^ as a function of N^°-' are shown 
for different colhding nuclei, In+In, S+S, Ne+Ne, 0+0 and C+C at Sjab^lSS AGeV. 

Fig. [16] presents the HSD scaled variances Wp^^ for C+C, 0+0, Ne+Ne, S+S, 
In+In, and Pb+Pb collisions at 158 AGeV as a function of Np'"-' . The fluctuations of 
N^p^^ are quite strong for mid-peripheral reactions with N^"-' = 20+30 and negligible for 
the most central collisions. A vanishing of u^p^^ = at Np^"-^ =A does not, however, show 
up in collisions of light nuclei (from C to S). Even for the maximal values of Np'"-^ = A 
the fluctuations Up^^ do not vanish and increase with decreasing atomic mass number 
A. For example in C+C collisions for Np'"-' = A = 12 the number of participants 
from the target still fluctuates and the scaled variance amounts to oop^^ = 0.25. Some 
combination of the system size Np and collision energy Eiab might move the chemical 
freeze-out point close the QCD critical point. One could then expect an increase of 
multiplicity fluctuations in comparison to their 'background values'. 

Why does one need central collisions of light and intermediate ions instead of 
studying peripheral Pb+Pb collisions for a search of the critical point? Fig. [T6] explains 
this issue. At fixed Np"-^ the average total number of participants, Np = Np'"-' +Np^'^ , is 
equal to (Np) = 2N''^°\ and, thus, it fluctuates as up = 0.5a;p^^. Then, for example, the 
value of N'p°-' = 30 corresponds to almost zero participant number fluctuations, up = 0, 
in S+S collisions while up becomes large and is close to 1 and 1.5 for In+In and Pb+Pb, 
respectively. Even if Np'"-' is fixed exactly, the sample of the peripheral collision events 
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in the heavy-ion case contains large fluctuations of the participant number: this would 
'mask' the critical point signals. As also seen in Fig. [16] [right), the picture becomes 
actually more complicated if the atomic mass number A is too small. In this case, the 
number of participants from a target starts to fluctuate significantly even for the largest 
and fixed value of Np'"-' =A. 

5.2. Multiplicity Fluctuations at Zero Impact Parameter 

The importance of a selection of the most central collisions for studies of hadron 
multiplicity fluctuations has been stressed in the previous sections. Due to its 
convenience in theoretical studies (e.g., in hydrodynamical models) one commonly uses 
the condition on impact parameter h for the selection of the 'most central' collisions in 
model calculations. However, the number of participant even at 6 = is not strictly 
fixed and fluctuates according to some distributions. It should be also stressed that the 
conditions h < bmax can not be fixed experimentally since the impact parameter itself 
can not be measured in a straightforward way. Actually, in experiments one accounts 
for the 1%, 2% etc. most central events selected by the measurement of spectators in the 
Veto calorimeter, which corresponds to the event class with the largest N'p"-' . As it will 
be demonstrated below the multiplicity fluctuations are very sensitive to the centrality 
selection criteria. In particular, the transport model results for 6 = and for 1% events 
with the largest N^"-' are rather different (see below). 



1.0 



0.9 



rii 0.8 



0.7 



0.6 



0.5 



HSD & UrQMD, b^O \ - 




HSD UrQMD 

— ■ ^ E=10 AGeV 

— • O- E=158AGeV 



40 80 120 160 200 
A, mass number 



HSD & UrQMD, b^O | , , 



HSD UrQMD 

d- E=10 AGeV 
-• O- E=158 AGeV 




40 80 120 160 
A, mass number 



200 



Figure 17. Left: Mean (Np), divided by the maximum number of participants 2A 
in events with 6 = for different nuclei at colhsion energies Eiab=^0 and 158 AGeV. 
Right: The scaled variance cjp in events with 6 = for different nuclei at collision 
energies Eiab=W and 158 AGeV. 



Let's start with the 6 = centrality selection criterium. First, we estimate the 
average number of participants, (Np), and the scaled variances of its fluctuations, up, 
in A + A collision events which satisfy the 6 = condition. The left panel in Fig. [T7I 
shows the ratio, (Np) /2A, in A + A collisions with 6 = for different nuclei at collision 
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energies Eiab = 10 and 158 AGeV. Both transport models (HSD and UrQMD) show 
a monotonous increase of (Np) /2A with colhsion energy for all nuclei in the energy 
range 10 158 AGeV (Fig. [T71 left). Correspondingly, the fluctuations of the number of 
participants up for all nuclei become smaller with increasing collision energy (Fig. [T71 
right). As seen from Fig.^7\{left) about 90% of all nucleons are participants for Pb+Pb 
collisions with b = 0. This number becomes essentially smaller, about 60-70%, for C+C 
collisions. One can therefore expect that participant number fluctuations at 6 = are 
small for heavy nuclei but strongly increase for light systems. This is demonstrated in 
Fig. HZ] (right): up is about O.l-i-0.2 in Pb+Pb and In+In but becomes much larger, 
0.5^0.7, in C+C collisions. 

One can conclude that the condition 6 = corresponds to 'most central' A + A 
collisions only for nuclei with large atomic mass number (In and Pb). In this case the 
average number of participants is close to its maximum value and its fluctuations are 
rather small. However, in the studies of event-by-event multiplicity fluctuations in the 
collisions of light nuclei (C and S) the criterium 6 = is far from selecting the 'most 
central' A + A collisions. 

Results of HSD and UrQMD transport model calculations for the scaled variance of 
negative, c<;_, positive, and all charged, Uch, hadrons are shown in Fig. [T8]at different 
collision energies, Eiab = 10, 20, 30, 40, 80, 158 AGeV, and for different colliding nuclei, 
C+C, S+S, In+In, Pb+Pb. The transport model results correspond to collision events 
for zero impact parameter, b = 0. To make the picture more complete, the transport 
model results for inelastic p+p collisions are shown too, for reference. Note that the 
proton spectators are not accounted for in the calculation of and Nch- Thus, proton 
spectators do not contribute to u+ and Uch- 

One sees a monotonic dependence of the multiplicity fluctuations on both Eiab and 
A: the scaled variances uj+, and Uch increase with Eiab and decrease with A. The 
results for p+p collisions are different from those for light ions. Note that within HSD 
a detailed comparison of the multiplicity fluctuations in N+N inelastic collisions and 
6 = heavy- ion collisions (Pb+Pb and Au+Au), including the energy dependence up to 
^/sNN = 200 GeV, has been presented in Sec. [Hand in Refs. [STtlFTj . Fig. [T8] corresponds 
to the full in acceptance. 

The combination of the multiplicity fluctuations in elementary N+N collisions 
and fluctuations of the number of nucleon participants explains the main features of 
hadron multiplicity fluctuations in A + A collisions. In particular, within the MIS 
f l20|) one obtains the dependence on collision energy and atomic mass number shown 
in Fig. [iHl The value of Ui in Eq. ( l20|) is the average number of i'th particles per 
participant, Ui = (Ni) /{Np), and cop equals the scaled variance for the number of 
nucleon participants. The N+N collisions deflne the fluctuations u* = ujf^ from a 
single source. 

In Fig. [19] the HSD results for Ui in A + A collisions at 6 = are compared to the 
Eq. ( 120|) . One concludes that the transport model results for the multiplicity fluctuations 
are in qualitative agreement with MIS (120|) . Both rii and u* increase strongly with 
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Figure 18. The results of HSD (left) and UrQMD (right) simulations for aj_ (top 
panel), a;+ (middle panel), and uJch (lower panel) in p+p and central C+C, S+S, 
In+In, Pb+Pb collisions at Eiab = 10, 20, 30, 40, 80, 158 AGeV. The condition 6 = 
is used here as a criterium for centrality selection. There are no cuts in acceptance. 

collision energy as seen from Fig. This explains, due to Eq. (l20l) . the monotonous 
increase with energy of the scaled variances Ui in A+ A collisions at 6 = seen in Fig. [181 
Note that up at b = decreases with collision energy as shown in Fig. [171 i^ight. This, 
however, does not compensate a strong increase of both rii and u*. The atomic mass 
number dependence of the scaled variances Ui in A + A collisions with 6 = follows 
from the A-dependence of cup. Fig. [T71 {right) demonstrates a strong increase of cjp for 
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Figure 19. The left panel illustrates the energy dependence of uj- in S+S collisions 
at 6 = in the full An acceptance, the right panel - the w_ dependence on atomic 
mass number at Eiab=SO AGeV. The HSD results are shown by the squares while the 
circles correspond to the MIS ((20)) . The stars show the first term, wl, in the r.h.s. of 
Eq. ((20)) - the scaled variance for negative hadrons in N+N collisions. The values of 
cj'H, rii, and cjp are calculated within HSD. 

light nuclei. This, due to ( 120)) . is transformed to the corresponding behavior of Ui seen 
in Fig. [m 

5.3. Centrality Trigger with the Number of Projectile Participants 

In this subsection the centrality selection procedure by fixing the number of projectile 
participants Np'"-' is concidered. This corresponds to the real situation of A+A collisions 
in fixed target experiments. As a first step one simulates in HSD and UrQMD the 
minimal bias events and calculate the event distribution over the number of participants 
Npart- Then, one selects 1% most central collisions which correspond to the largest values 
of Np'"-' . In such a sample of A+A collisions events with largest Np'"-' from different 
impact parameters can contribute. After that one calculates the values of in these 
samples. Note, that even for a fixed number of Np'"-^ the number of target participants 
Np^^ fluctuates. Thus, the total number of participants Np = Np^^ + Np'°-' fluctuates, 
too. In our 1% sample, both Np'^^ and Np'"-' fluctuate. Besides there are correlations 
between iVp^^ and Np'"-' . 

Fig. 120) shows the ratio {Np)/2A and the scaled variance, up, for 1% most central 
collisions selected by the largest values of N^°-' . These results are compared with those 
for the 6 = centrality selection. For heavy nuclei, like In and Pb, one finds no essential 
differences between these two criteria of centrality selection. However, the 1% centrality 
trigger defined by the largest values of N^°-' looks much more rigid for light ions (S and 
C). In this case the ratio {Np)/2A is larger, and is essentially smaller than for the 
criterion 6 = 0. As a result the 1% centrality trigger by the largest values of Np'"-' leads 
to a rather weak A-dependence of up. 




Figure 20. The HSD {left) and UrQMD (right) results for the ratfo {Np)/2A (the 
upper panel) and the scaled variance of the participant number fluctuations, ujp (the 
lower panel), for the 1% most central collisions selected by the largest values of N^°-' 
(full symbols), for different nuclei at collision energies i?;a6=10 and 158 AGeV. The 
open symbols present the results of Fig. [T7] [right) for h — Q. 



Some comments are appropriate at this point. Let's define the centrahty c{N) 
as a percentage of events with a multiplicity larger than (this can be the number 
of produced hadrons, number of participants, etc.). It was argued in Ref. [72] that a 
selection of c{N) of most central A + A collisions is equivalent to restricting the impact 
parameter, h < b{N), with, 

b(N) = c(iV) , (29) 

V 71 

where aimi is the total inelastic A+A cross section. Thus, the centrality criterion by 
the multiplicity A^ is equivalent to the geometrical criterion by the impact parameter b. 
Moreover, the result (1291) does not depend on the specific observable A^ used to define 
the c-percentage of most central A+A collisions. The result f l29|) should be valid for any 
observable A^ which is a monotonic function of b. Therefore, the relation f l29|) reduces 
any centrality selection to the geometrical one. This result was obtained in Ref. [72] by 
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neglecting the fluctuations of multiplicity at a given value of b. This is valid if c is not 
too small and the colliding nuclei are not too light. In the sample of A+A events with 
1% of largest N^°\ the relation ( 129|) can not be applied for S+S and C+C collisions. 
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Figure 21. The dependence of uj- on atomic mass number at i?iab=80 AGeV for 
the HSD {left) and UrQMD (Ze/f) simulations. The squares correspond to & = 0, and 
circles to 1% largest N^°-' . 



Fig. EH shows a comparison of the A-dependence of in the transport models for 
two different samples of the collision events: for 6 = and for the 1% of events with 
largest N^"-' values. One can see that the multiplicity fluctuations are rather different in 
these two samples. Moreover, these differences are in the opposite directions for heavy 
nuclei and for light nuclei. For light nuclei, a;_ is essentially smaller in the 1% sample 
with largest N^°-' values, whereas for heavy nuclei the smaller fluctuations correspond 
to 6 = events. Note that in the 1% sample with largest N'p°-' values the A-dependence 
of multiplicity fluctuations becomes much weaker. In this case a strong increase of the 
multiplicity fluctuations for light nuclei, seen for 6 = 0, disappears. 

For the 1% most central A + A collision events - selected by the largest values of 
j\fpro] _ HSD multiplicity fluctuations are shown in Fig. |22] and Fig. |23l For light 
nuclei (S and C) the multiplicity fluctuations in the samples of 1% most central collisions 
are smaller than in the 6 = selection and the atomic mass number dependencies become 
less pronounced (compare Fig. [22] and Fig. ITSll . This is because the participant number 
fluctuations up have now essentially smaller A-dependence, as seen in Fig. |20l 

Fig. HE] shows that both HSD and UrQMD predict a monotonic dependence of the 
charge particle multiplicity with energy. So, the hadronic 'background' for the NA61 
experiments is expected to be a smooth monotonic function of beam energy. 

The event-by-event observables show a higher sensitivity to the initial nucleon 
density distribution than the standard single particle observables [73]. A sensitivity 
of the A-dependence of cuch to these details of the models indicates a necessity for 
further studies of the initializations of the nuclei in transport model approaches. This 
becomes important for the theoretical interpretation of future experimental data on 
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Figure 22. The HSD results for a;_ (upper panel), (middle panel), and ujch (lower 
panel) in A + A and p+p collisions for the full Att acceptance in 3D (left) and 2D 
(right.) projection. The 1% most central C+C, S+S, In+In, and Pb+Pb collisions are 
selected by choosing the largest values of N^°-' at different collision energies i?;ab=10, 
20, 30, 40, 80, 158 AGeV. The errorbars indicate the estimated uncertainties in the 
model calculations. The HSD results from inelastic p+p collisions are the same as in 
Fig.IH 



event-by-event fluctuations. 

Note that the MIS and Eq. (I20l) work for the multiphcity fluctuations simulated 
by the transport models in full 47r acceptance but not for the acceptance in a speciflc 
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Figure 23. The same as in Fig. [521 but for final hadrons accepted in the projectile 
hemisphere, y > (left), and in the target hemisphere, y < (right). 

rapidity region. The results for inelastic p+p collisions are identical in the projectile and 
target hemispheres. This is not the case in the sample of 1% most central A+A collisions 
selected by N^°-' . The total number of nucleons participating in A + A collisions 
fluctuates. These fluctuations are not symmetric in forward-backward hemispheres: 
in the selected 1% sample the number of target participants N^p^^ fluctuates essentially 
stronger than that of N^°-' . The HSD results in Fig. [23] clearly demonstrate larger 
values for all scaled variances, W-, uj+, and Uch, for ?/ < acceptance than those for 
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y > one. This is due to stronger target participant fluctuations, Up^^ > Wp"-'. 

6. Fluctuations and Correlations in Au+Au Collisions at RHIC 

In this Section we discussed some features of the particle number fluctuations ^74j and 
forward-backward correlations [75] in Au+Au collisions at ^/sNN = 200 GeV. 

6.1. Multiplicity Fluctuations in Au+Au Collisions at RHIC 

The charged multiplicity fluctuations in Au+Au collisions at ^/sNN = 200 GeV have 
been measured recently by the PHENIX Collaboration [7^177] . The centrality selection 
is an important aspect of fluctuation studies in A+A collisions. As it was discussed in the 
previous sections the samples of collisions with a fixed number of projectile participants 
Np'"-' can be selected to minimize the participant number fluctuations at the fixed target 
experiments. This selection is possible due to a measurement of the number of nucleon 
spectators from the projectile, A|'^°-', in each individual collision by a calorimeter which 
covers the projectile fragmentation domain. In a collider type experiments another 
centrality trigger should be used. For the multiplicity fluctuations data the PHENIX 
Collaboration uses two kinds of detectors which define the centrality of Au+Au collision: 
Beam-Beam Counters (BBC) and Zero Degree Calorimeters (ZDC). The BBC measure 
the charged particle multiplicity in the pseudorapidity range 3.0 < \ri\ < 3.9, and the 
ZDC - the number of neutrons with \ri\ > 6.0 [76l[77]. These neutrons are part of 
the nucleon spectators. Due to technical reasons the neutron spectators can be only 
detected by the ZDC (not protons and nuclear fragments), but in both hemispheres. 
The BBC distribution will be used in the HSD calculations to divide Au+Au collision 
events into 5% centrahty samples. The HSD does not specify different spectator groups - 
neutrons, protons, and nuclear fragments. Thus, one can not use the ZDC information. 
In Fig. [21] (left) the HSD results for the BBC distribution and centrality classes in 
Au+Au collisions at ^^snn =200 GeV are shown. One finds a good agreement between 
the HSD shape of the BBC distribution and the PHENIX data [7S1[77|. Note, however, 
that the HSD (Np) numbers are not exactly equal to the PHENIX values. It is also not 
obvious that different definitions for the 5% centrality classes give the same values of 
the scaled variance up for the participant number fluctuations. 

Defining the centrality selection via the HSD transport model (which is similar 
to the BBC in the PHENIX experiment) we calculate the mean number of nucleon 
participants, {Np), and the scaled variance of its fluctuations, oup, in each 5% centrality 
sample. The results are shown in Fig. |2ll right. The Fig. ^5\{left) shows the HSD results 
for the mean number of charged hadrons per nucleon participant, rij = (Ni) / (Np) , 
where the index i stands for "— ", "+", and "ch" final hadrons. Note that the centrality 
dependence of Ui is opposite to that of cop. 

The PHENIX detector accepts charged particles in a small region of the phase space 
with pseudorapidity \ri\ < 0.26, azimuthal angle (p < 245°, and the pt range from 0.2 to 
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Figure 24. HSD model results for Au+Au collisions at y^s^N — 200 GeV. Left: 
Centrality classes defined via the BBC distribution. Right: The average number of 
participants, (-/Vp), and the scaled variance of the participant number fluctuations, 
ujp, calculated for the 5% BBC centrality classes. 




Figure 25. HSD results for different BBC centrality classes in Au+Au collisions 
at ^ysNN — 200 GeV. Left: The mean number of charged hadrons per participant, 
n, = {N,)/{Np). Right: The fraction of accepted particles, qt = {N"-"") / (N,) . 

2.0 GeV/c [ISIIIT]. The fraction of the accepted particles Qi = {Nf''^)/{Ni) calculated 
within the HSD model is shown in Fig. [251 right. According to the HSD results only 
about 3% of charged particles are accepted by the mid-rapidity PHENIX detector. 

To estimate the role of the participant number event-by-event fluctuations the 
MIS relation ( 120|) has been used. It has been assumed that N-l-N collisions define the 
fluctuations uj* from a single source. To calculate the fluctuations u'^'^^ in the PHENIX 
acceptance the acceptance scaling formula (128!) is used. One finds, 

' = I - q, +q,uj* + Qi n, cop . (30) 

The HSD results for up (Fig. El right), rii (Fig. ESI left), qi (Fig. EHl right), 
together with the HSD nucleon-nucleon values, cul = 3.0, u*^ = 2.7, and cu*^ = 5.7 at 
^/sNN = 200 GeV, define completely the results for u"''^'^ according to ( l30l) . We find a 
surprisingly good agreement of the results given by (l30l) with the PHENIX data shown 
in Fig. EH The centrality dependence of u"''^'^ stems from the product, qtUiUp, in the 
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Figure 26. The scaled variance of charged particle fluctuations in Au+Au collisions 
at ,/sWn = 200 GeV with the PHENIX acceptance. The circles are the PHENIX 
data [76l[77] while the open points (connected by the solid line) correspond to pO|) 
with the HSD results for wp, rii, and qi. 

last term of the r.h.s. of f pOj) . Note that up decreases with (Np), whereas both rii 
and Qi increase. It may lead to a nontrivial centrality dependence with a maximum at 
intermediate values of (Np) as seen in the PHENIX data. 

6.2. Forward-Backward Correlations 

Correlations of particles between different regions of rapidity have for a long time been 
considered to be a signature of new physics. A shortening in the correlation length 
in rapidity has been thought to signal a transition to a quark-gluon plasma [78|I79]. 
Conversely, the appearance of long-range correlations has been associated with the onset 
of the percolation limit, also linked to the QCD phase transition [SniEl]- Recently, the 
correlations across a large distance in rapidity have also been suggested to arise from a 
color glass condensate |H2l[H3]. The observation of such correlations in A+A collisions 
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at RHIC energies by the STAR Collaboration (HUES] has therefore elicited a lot of 
theoretical interest. 

The purpose of this Section is to identify some baseline contributions to the 
experimentally observed correlations, contributions that do not depend on new physics 
[75] . Two models that incorporate event-by-event fluctuations in initial conditions have 
been used to illustrate the effect of these contributions: the HSD transport model and 
a 'toy' wounded nucleon model. Based on such a comparison one can argue that a 
study of the dependence of correlations on the centrality bin definition as well as the 
bin size may distinguish between 'trivial' correlations and correlations arising from 'new 
physics'. 

The statistical properties of a particular sample of events can be characterized by 
a set of moments or cumulants of some observable. These properties depend upon a 
set of criteria which are used to select this sample. Applied to the context of A+A 
collisions this translates to the construction of centrality bins of collision events from 
minimum-bias data. The charged hadron multiplicities Na and Nb will be considered 
in two symmetric intervals Arj of pseudo-rapidity. After construction of the centrality 
bins, one can calculate the moments of a resulting distribution P^^''p{Na-, Nb', Arj): 

{Nl ■ AT^j)^-^ = ^ ivj^ P,^«-(Ar^, Nb; Ar/) . (31) 

Na,Nb 

In fl3T|) the subscript c denotes a particular centrality bin, while the superscript rjgap 
denotes the separation of two symmetric intervals Ar] in pseudo-rapidity space where 
particle multiplicities Na and Nb are measured. The correlation coefficient is defined 
by 

and measures how strongly multiplicities Na and Nb - in a given centrality bin c for 
pseudo-rapidity separation rjgap - are correlated. In (1521) . AA^ = N — (A^)^""'' and 
((AA'^)^)^'"''' = {{ANb)'^)^^"^ for symmetric intervals. 

The recent preliminary data on the forward-backward correlation coefficient (l32l) of 
charged particles by the STAR Collaboration [5^155] exhibit two striking features: a) an 
approximate independence on the width of the pseudo-rapidity gap rjgap , b) a strong 
increase of p with centrality. 

6.3. Glauber Monte Carlo Model 

The PHOBOS Glauber Monte Carlo code [86] coupled to a 'toy' wounded nucleon model 
is used here, referred to as CMC. The aim of this model is to emphasize two crucial 
aspects: 1) an averaging over different system sizes within one centrality bin introduces 
correlations; 2) the strength of these correlations depend on the criteria used for the 
centrality definition and on the size of the centrality bins. 

Employing the Glauber code the distribution of the number of participating 
nucleons, Np, is modelled in each A+A collision for given impact parameter b (cf. 
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Fig. |271 left). This is done for Au+Au with standard Wood-Saxon profile and the N+N 
cross section of a^^ = 42 mb. The 'event' construction proceeds then in a two-step 
process. Firstly, the total number of charged particles is randomly generated: 

Np 

Nch = Y.<h. (33) 

i=l 

where the number of charged particles n^^ per participating nucleon are generated 
by independently sampling a Poisson distribution with given mean value rich = 10. 
Secondly, these charged particles are randomly distributed according to a Gaussian in 
pseudo-rapidity space: 

where cr^ = 3 defines the width of the pseudo-rapidity distribution. Hence, in each single 
event there are no correlations between the momenta of any two particles. Note that 
numerical values of rich and are fixed in a way to have a rough correspondence with 
the data on charged particle production at yJSNN = 200 GeV. 




Figure 27. Left: The histogram shows the distribution of events with a fixed number 
of participating nucleons Np and fixed impact parameter b in Au+Au collisions at 
y^SATjv — 200 GeV. Right: The scaled variance cjp of the distribution of participating 
nucleons in 10% bins as defined via b, Np, and N^l^^ . 

Fig. 127] (left) shows the CMC event distribution in the (6, A'^p)-plane. For each 
of these events we randomly generate the number of charged particles Nch and their 
?7-distribution according to (133|) and (IMj) . respectively. The construction of centrality 
classes can now be done in several ways. The following criteria are chosen: via impact 
parameter b, via the number of participating (wounded) nucleons Np, and via the 
charged particle multiplicity N^j^^ in the midrapidity window \ri\ < 1. 

In the case of choosing the number of participating nucleons Np for centrality 
definition, one takes vertical cuts in Fig. |27] (/e/f), while choosing the impact parameter 
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b, one takes horizontal cuts. Hence, depending on the centrahty definition, one may 
assign a particular event (characterized by Np and b) to two different centrality bins. 
Fig. ET] (right) shows the resulting scaled variance cop, 

{{ANp)\ 



OOp = 



{Np)c 



(35) 



of the underlying distribution of the number of participating nucleons Np in each 
centrality bin. Using the centrality selection via impact parameter b, which is only 
the theoretically available trigger, one generally obtains a rather wide distribution of 
participating nucleons in each bin. The lines for centrality selections via N^^^^ and via 
Np are similar due to the event construction by (133|) and ([3i 
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Figure 28. The forward-backward correlation coefficient p for 10% centrality classes 
defined via Np (left), via the impact parameter b (middle), and via the multiplicity in 
the central rapidity region N^^J (right). 



The sensitivity of the forward-backward correlation signal as a function of the 
separation rjg^p of two narrow intervals (At] = 0.2) on the centrality definition is 
investigated now. This is done for the 10% centrality defined via Np, via b, and via 
N^^. The results are shown in Fig. [281 In the GMC one can identify the number of 
participating nucleons Np with the system size, and as the measure for system size 
fiuctuations. Having a large system - as measured by Np - implies a large number 
of charged particles Nch- In GMC they are distributed independently in pseudo- 
rapidity space. Conversely, an event with small Np contains only few charged particles. 
By grouping the collision events into 10% centrality bins one finds rather large Np- 
fluctuations in one specific bin. The averaging over different states in the centrality 
bin introduces correlations between any two regions of pseudo-rapidity. Small systems 
will have few particles 'on the left' and few particles 'on the right' with respect to 
midrapidity. Large systems will have many particles 'on the left' and many particles 
'on the right'. But this just means a non-zero forward-backward correlation. From the 
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definition (132 p one finds a positive correlation coefficient p due to averaging over system 
sizes. 

Note that centrality selections via Np and via Nch give essentially the same results 
for p in the GMC (cf. left and right panels of Fig. [28|) . Using the impact parameter b 
for the centrality definition generates centrality bins with almost constant p as seen in 
Fig. {middle). This is due to a rather fiat dependence of cop on the centrality defined 
via b as shown in Fig. [27| {right). In the GMC model the apparent ordering of p values 
with respect to centrality bins originates from the width of the underlying distribution 
in the number of wounded nucleons in each bin, i.e. from the values of up . 

The measured and apparently strong forward-backward correlations can be 
accounted for by a 'toy' model such as the GMC, provided it produces particles over 
the whole rapidity range and includes strong enough event-by-event fluctuations of Np. 
The next section will show that an introduction of dynamics and hadron re-interactions 
within HSD does not alter these conclusions significantly. 

6.4- HSD Transport Model Simulations 

A physically more reasonable scenario which, however, also does not include any 'new 
physics' (such as color glass condensate, quark-gluon plasma, etc.) can be obtained in 
the HSD transport approach. As before within the GMC, the HSD events are generated 
according to a uniform distribution, Nev{b) ~ b. The resulting distribution of events in 
the {Np, 6)-plane is similar to the GMC result depicted in Fig. [27] {left). Note, that the 
peripheral part of the distribution determines also the centrality binning and the real bin 
widths. This is crucial for most central collisions where the number of events is small. 
Slight uncertainties in the peripheral "tail" of the distribution leads to large errors in 
the sizes of most central bins and hence to large changes in results for fluctuations and 
correlations. 

In contrast to the STAR data, the charged particle reference multiplicity iV^^-^ in the 
same pseudo-rapidity range |?7| < 1 for all values of rjgap is used in the HSD simulations. 
This procedure introduces a systematic bias, since the pseudo-rapidity regions for the 
measured multiplicity in a small Arj window (signal) and for the reference multiplicity 
partially overlap. This bias, however, is small and does not affect any of the conclusions. 

Fig. [291 shows the scaled variance of the underlying Np distribution for 10% {left) 
and 2% {right) centrality bins defined via different centrality triggers within HSD. The 
results for 10% bins can be compared with the scaled variance up in the GMC model 
in Fig. [27] {right). Fluctuations of the number of participants, as well as their average 
values, are similar in both HSD and GMC models when the centrality bins are defined 
via Np. These quantities are completely defined by the Np distribution, which is similar 
in both models. Binning via the impact parameter b in HSD, as well as in GMC, gives 
decreasing fluctuations in the participant number with increasing collision centrality. 
The results for 10% bins defined via the reference multiplicity are rather different in the 
GMC and HSD models. In GMC the charged multiplicity distribution is implemented 
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Figure 29. The HSD results for the fluctuatfons uipart as a function of the mean value 
(Np) of the participating nucleons within bins as defined via 6, Np, and N^^ . The 
left panel corresponds to a 10% and the right to a 2% bin width. 

according to ( 133|) and (IMj) . Hence, the results obtained by binning via the reference 
muhiphcity follow the line obtained by binning via Np. In contrast to the GMC, in the 
HSD simulations the average number of charged particles rich per participating nucleon 
is not a constant, but increases with Np. Additionally, the shape of rapidity distribution 
is also different in different centrality bins. These two effects lead to different values of 
cup in the centrality bins defined via N^^^ in the GMC and HSD models. An interesting 
feature seen in Figs. ET] and [30] is that up is large for 10% most central bins. The 
behavior for 2% centrality bins is rather different: up decreases with centrality as seen 
from Fig. [29] right. This becomes similar to the dependence of up on Np'"-^ considered 
in Sec. [3] for the fixed target experiments. 

One comment is appropriate here. It was argued in Ref. [72] that any centrality 
selection in A+A collisions is equivalent to the geometrical one via impact parameter 
b. This was already discussed in Sec. [1] Different centrality selection criterions give 
indeed the same average values of physical observables. However, they may lead to 
rather different fluctuations of these observables in the corresponding centrality bins, 
cf. equal values of (Np) and different values of for different centrality selections 
presented in Fig. [29] When considering smaller centrality bins (2% in Fig. [29] right) the 
fluctuations in the participant number become smaller but more strongly dependent on 
the definition of the binning. 

Fig. [30l summarizes the dependence of forward-backward correlation coefficient p as 
a function of rjgap on the bin size and centrality definition within the HSD model. The 
dependence of p on rjgap is almost fiat, reflecting a boost-invariant distribution of particles 
created by string breaking in the HSD. The right top panel of Fig. [30] demonstrates also a 
comparison of the HSD results with the STAR data [8H[85] . One observes that the HSD 
results exceed systematically the STAR data. However, the main qualitative features 
of the STAR data - an approximate independence of the width of the pseudo-rapidity 
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Figure 30. The HSD results for the forward-backward correlation coefficient p for 10% 
(top) and 2% (bottom) centrality classes defined via Np (left), via impact parameter b 
(center), and via the reference multiplicity N^l^'^ (right). The symbols in the top right 
panel present the STAR data in Au+Au collisions at y^SNN = 200 GeV [84l|85] . 



gap r]gap and a strong increase of p with centrality - are fully reproduced by the HSD 
simulations. 

The correlation coefficient p largely follows the trend of the participant number 
fluctuations as a function of centrality. The actual results, however, strongly depend 
on the way of defining the centrality bins. For instance, choosing smaller centrality bins 
leads to weaker forward-backward correlations, a less pronounced centrality dependence, 
and a stronger dependence on the bin definition. The physical origin for this is 
demonstrated in Fig. [31] As the bin size becomes comparable to the width of the 
correlation band between Np and N^^'^ , the systematic deviations of different centrality 
selections become dominant: the same centrality bins defined by Np and by iVj^'^ contain 
different events and may give rather different values for the forward-backward correlation 
coefficient p. 

It should be underlined that these properties are specific to the geometric nature 
of the correlations analyzed here. If the observed fiuctuations are of dynamical origin 
(for example, arising from the quantum fiuctuations of coherent fields created in the 
first fm/c of the system's lifetime as in Refs. [82l[83]), there are no evident reasons why 
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Figure 31. The histogram shows the distribution of HSD events with fixed number of 
participating nucleons Np and fixed reference charge particle multiplicity N^jJ ■ The 
same centrality class (20-22% as an example) defined in various ways contains different 
events. 

they should strongly depend on the centrality bin definitions and bin sizes. Thus, the 
experimental analysis for different bin sizes and centrality definitions - as performed 
here - may serve as a diagnostic tool for an origin of the observed correlations. A strong 
specific dependence of the correlations on bin size and centrality definition would signify 
their geometrical origin. 

7. Fluctuations of Conserved Charges in Nucleus-Nucleus Collisions 

The aim of the present Section is to study event-by-event fiuctuations of the conserved 
charges - net baryon number and electric charge - in Pb-Pb collisions at 158 AGeV 
energies within the HSD transport approach [87] . 

7.1. Net Baryon Number Fluctuations 

Let's start with a quantitative discussion by first considering the fiuctuations of the 
net baryon number in different regions of the participant domain in collisions of two 
identical nuclei. These fluctuations are most closely related to the fluctuations of the 
number of participant nucleons because of baryon number conservation. 

The HSD results for lob in Pb+Pb at 158 AGeV are presented in Fig. [32l In 
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Figure 32. The HSD simulations for Pb+Pb collisions at 158 AGeV for fixed values 
of Np'^"-'. Left: The baryon number fluctuations in full acceptance, w^, in projectile 
hemisphere, ujg (lower curve), and in target hemisphere, ujg (upper curve). Right: 
The scaled variances of the baryon number fluctuations in different rapidity intervals. 



each event the nucleon spectators when counting the number of baryons are subtracted. 
Otherwise there would be no fluctuations of the net baryon number in the full acceptance 
and ub equals to zero. The net baryon number in the full phase space, B = Nb — N-^, 
equals to the total number of participants Np = N^p^^ + N^"-' . At fixed Np^"-' the 
Np number fluctuates due to fluctuations of Np^^ . These fluctuations correspond to 
an average value, {N^p"^^) ~ N^"-' , and a scaled variance, oo^p^^ (see Fig. H]). Thus, for 
the net baryon number fluctuations in the full phase space one finds, ub — 0.5up^^. A 
factor 1/2 appears because only a half of the total number of participants fluctuates. 

One introduces Up and Up, where the superscripts p and t mark quantities measured 
in the projectile and target momentum hemispheres, respectively. Fig. [32] demonstrates 
that Ub > both in the whole projectile-target hemispheres and in the symmetric 
rapidity intervals. On the other hand one observes that Ub ~ ojb iii most central 
collisions. This is because the fluctuations of the target participants become negligible 
in this case, i.e. Wp^^ = (Fig. HI right). As a consequence the fluctuations of any 
observable in the symmetric rapidity intervals become identical in most central collisions. 
Note also that transparency-mixing effects are different at different rapidities. From 
Fig. [32] [right) it follows that Up in the target rapidity interval [—2, —1] is much larger 
than Up in the symmetric projectile rapidity interval [1, 2]. This fact reveals the strong 
transparency effects. On the other hand, the behavior is different in symmetric rapidity 
intervals near the midrapidity. From Fig. [32] (right) one observes that Up in the target 
rapidity interval [—1, 0] is already much closer to Up in the symmetric projectile rapidity 
interval [0, 1]. This gives a rough estimate of the width. Ay ^ 1, for the region in rapidity 
space where projectile and target nucleons communicate to each others. 
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7.2. Net Electric Charge Fluctuations 



It is difficult to study experimentally the baryon number fluctuations since an 
identification of neutrons in a large acceptance in a single event is difficult. In this 
subsection we consider the HSD results for the net electric charge, Q, fluctuations. As 
Q = 0.4B in the initial heavy nuclei one can naively expect that Q fluctuations are quite 
similar to B fluctuations. However, there is a principal difference between Q and B in 
relativistic A+A collisions. Fig. |33] demonstrates the rapidity distributions of the net 
baryon number, B = Nb — N-^ i^^ft), and total number of baryons, Nb + N-^ (right), for 
different centralities in Pb+Pb collisions at 158 AGeV. One observes that both quantities 
are very close to each other; the y-dependence and absolute values are very close for B 
and Nb + N-^ distributions. This is, of course, because the number of antibaryons is 
rather small, N-^ ^ Nb- 
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Figure 33. The HSD rapidity distributions in Pb+Pb collisions at 158 AGeV for the 
net baryon number, B = Nb — N-^ (left), and total number of baryons, Nb + N^ 
(right), at different N^°-' and in the minimum bias (m.b.) sample. 

Fig. IMl shows the same as Fig. [33] but for the electric charge Q = N+ — N_ (left), 
and total number of charged particles, Nch = N^ + N_ (right). The y-dependence of 
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Figure 34. The same as in Fig. [33] but for the electric charge Q — N+ — N_ (left), 
and total number of charged particles, Nch = N+ + N- (right). 



dQ/dy and dNch/dy is quite different. Besides, the absolute values of Nch are about 
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10 times larger than those of Q. This imphes that Q ^ ~ A^_. In the previous 
subsection the scaled variance ub to quantify the measure of the net baryon fluctuations 
has been used. It appears to be a useful variable as ojb is straightforwardly connected 
to u!p^^ and due to the relatively small number of antibaryons. Fig. [M] tells that uq is 
a bad measure of the electric charge fluctuations in high energy A+A collisions. One 
observes that uq = Var{Q)/ (Q) is much larger than 1 simply due to the small value of 
(Q) in a comparison with iV+ and N^. If the A+A collision energy increases, it follows, 
(A^_^) = (A^„) increases with collision energy, whereas (Q) = OAB remains constant. 
This leads to very large values of uq — The same would happen with ub, too, but at 
much larger collision energies. A useful measure of the net electric charge fluctuations 
is the quantity (see, e.g., [88]): 

Var{Q) 



X, 



Q 



(36) 



A value of Xq can be easily calculated for the Boltzmann ideal gas in the grand canonical 
ensemble. In this case the number of negative and positive particles fluctuates according 
to the Poisson distribution (i.e. = uj+ = 1), and the correlation between and 
N_ are absent (i.e. (N+N^) = (iV+)(iV_)), so that Xq = 1. On the other hand, the 
canonical ensemble formulation (i.e. when Q = const fixed exactly for all microscopic 
states of the system) leads to Xq = 0. Fig. [351 shows the results of the HSD simulations 
for the full acceptance, for the projectile and target hemispheres [left), and also for 
symmetric rapidity intervals in the c.m.s. (right). 
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Figure 35. Left: The HSD simulations in Pb+Pb collisions at 158 AGeV for Xq at 
different values of N^°-' in the full acceptance (lower curve) , for the projectile (middle 
curve) and target (upper curve) hemispheres. Right: The same, but for the fixed 
rapidity intervals in the c.m.s. 

The Q fluctuation in the full acceptance is due to N^p^^ fluctuations. As Q = OAB 
in colliding (heavy) nuclei, one may expect Var{Q) = 0.16 Var{B). In addition, 
(A^^^) ^ 4:{Np) at 158 AGeV, so that one estimates Xq = 0.04 cob for the fluctuations 
in the full phase space. The actual values of Xq presented in Fig. |35] (left) are about 3 
times larger. This is because of Q fluctuations due to different event-by-event values of 
proton and neutron participants even in a sample with fixed values of Np'°-' and N^p^^ . 
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From Fig. (right) one sees only a tiny difference between the Xq values in 
the symmetric rapidity intervals in the projectile and target hemispheres, and slightly 
stronger effects for the whole projectile and target hemispheres (Fig. [35], left). In fact, 
the fluctuations of and are very different in the projectile and target hemispheres, 
and the scaled variances and have a very strong A^p''°-'-dependence. This is shown 
on top- left and middle- left panels of Fig. |8l 

The Xq can be presented in two equivalent forms 

(iV+) , (iV_) ,A(iV+,iV_) 



X, 



Q 



2 UJ+ 



2 u 



(N-) 



- 2 



^ch 



(37) 



The Eq. ( l37ll is valid for any region of the phase space: full phase space, projectile or 
target hemisphere, etc. As seen from Fig. [HI both and are large and strongly 
A^p^°-^-dependent. This is not seen in Xq because of strong correlations between A^^ and 
N^, i.e. the term 2 A{N^, N_) / {Nch) compensates oj+ and uj- terms in (1371) . 

Fig. [36] shows a comparison of the HSD results for Xq with NA49 data in Pb+Pb 
collisions at 158 AGeV for the forward rapidity interval 1.1 < y < 2.6 inside the 
projectile hemisphere with additional p^-filter imposed. 
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Figure 36. The HSD results for Xq for Pb+Pb collisfons at 158 AGeV for the forward 
rapidity interval 1.1 < y < 2.6 inside the projectile hemisphere. The solid dots are the 
estimates obtained from Eq. (|37l) using the NA49 experimental data [46l|89l|9^ (the 
error bars are not indicated here). For illustration, the HSD results in the symmetric 
backward rapidity interval —2.6 < y < —1.1 (target hemisphere) are also presented. 

As an illustration, the HSD results in the symmetric backward rapidity interval 
—2.6 < y < —1.1 (target hemisphere) are also included. One observes no difference 
between the Xq results for the NA49 acceptance in the projectile and target hemispheres. 
The HSD values for uj+, and cUch are rather different in the projectile and target 
hemispheres for the NA49 acceptance (see Fig. [8]). This is not seen in Fig. [36] for Xq. 
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As explained above a cancellation between and lOch terms take place in fl37p . In 

fact, NA49 did not perform the Xq measurements. The Xg-data (solid dots) presented 
in Fig. [36] are obtained from Eq. (1371) using the NA49 data for 0;+, c<;_, and ujch as 
well as (A^+), (A^-), and {Nch) |l6l|89l[90]. Such a procedure leads, however, to very 
large errors for Xq (not indicated in Fig. 1361) which excludes any conclusion about the 
(dis) agreement of HSD results with NA49 data. 




Figure 37. The HSD results for electric charge fluctuations in 2% most central Pb+Pb 
collisions at 158 AGeV in the symmetric rapidity interval, AF = [— y, y] as a function 
of Ay — AF/2 in the c.m.s. A left panel shows the behavior of Xq, and a right one 
demonstrates separately and uj ch- 

it is instructive to consider the electric charge fluctuations in the symmetric rapidity 
interval [—y-,y\ in the c.m.s. for the most central Pb+Pb events. The sample of most 
central events is chosen by restricting the impact parameter to 6 < 2 fm. It gives about 
2% most central Pb+Pb collisions from the whole minimum bias sample. Fig. [371 shows 
the HSD results for electric charge fluctuations in 2% most central Pb+Pb collisions 
for the symmetric rapidity interval AY = [—y,y] in the c.m.s. as the function of 
Ay = AY/2. 

For Ay — )■ one finds Xq — )■ 1. This can be understood as follows: For AY — )■ 
the fluctuations of negatively, positively and all charged particles behave as for the 
Poisson distribution: u+ = co- = Uch — 1- Then from Eq. (1371) it follows that Xq = 1, 
too. From Fig. |371 (right) one observes that u^, and iOch all increase with increasing 
interval AY. However, Xq decreases with AY and - because of global Q conservation 
- it goes approximately to zero when all final particles are accepted. 

8. Fluctuations of Ratios in Nucleus-Nucleus Collisions 

The measurement of the fluctuations in the kaon to pion ratio by the NA49 Collaboration 
|91] was the first event-by-event measurement in nucleus-nucleus collisions. It was 
suggested that this ratio might allow to distinguish events with enhanced strangeness 
production attributed to the QGP phase. Nowadays, the excitation function for this 
observable is available in a wide range of energies: from the NA49 collaboration [92] 
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in Pb+Pb collisions at the CERN SPS and from the STAR collaboration in 
Au+Au collisions at RHIC. First statistical model estimates of the K/tt fluctuations 
have been reported in Refs. [951196] . The results from the transport model UrQMD have 
been presented in Ref. [97] for the SPS energy range. 

This Section presents the results of a systematic study of K/ti, K/p and p/vr ratio 
fluctuations [K = + K~ , vr = 7r+ + vr", and p means p + p) based on the HSD 
transport model. These results will be copmared to statistical model results in different 
ensembles [98l[99] and available experimental data. 



8.1. Notations and Approximations 

The fluctuations of the ratio Rab = Na/Nb will be characterized by [95] 
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(38) 



Using the expansion, 

Na ^ 
Nb 



{Na) + ANa _ {Na) + ANa 
{Nb) + ANb 

ANb 



AN, 

{Nb) ' \{Nb) 




X 



(39) 



one finds to second order in ANa/{Na) and ANb/{Nb) the average value and 
fiuctuations of the A to B ratio: 
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{Na) {Nb) [{Na){Nb) 
If species A and B fiuctuate independently according to the Poisson distributions 
(this takes place, for example, in the GCE for an ideal Boltzmann gas) one finds, 
ua = <^B = ^ and pab = 0. The Eq. (14T|) then reads 

1 1 



In a thermal gas, the average multiplicities are proportional to the system volume V. 
The Eq. (H2|) demonstrates then a simple dependence a' cx 1/V on the system volume. 

A few examples concerning to Eq. (HTl) are appropriate here. When {Nb) ^ {Na), 
e.g., A = + K~ and B = + , the a' (jUj) is dominated by the less abundant 
particles and the resonances decaying into it. When {Na) = {Nb), e.g., A = and 
B = 71^, the correlation term in Eq. f HTl) may become especially important. A resonance 
decaying always into a vr+Tr^-pair does not contribute to a' ( HTj) . but contributes to vr"*" 
and 7T~ average multiplicities. This leads [96] to a suppression of a' (HTl) in comparison 
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to its value given by Eq. f l42p . For example, if all tt"*" and vr" particles come by pairs 
from decay of resonances, one finds the correlation coefficient pj^+j^- = 1 in Eq. f HTj) . 
and thus = 0. In this case, the numbers of vr"*" and vr" fluctuate as the number of 
resonances, but the ratio vt+Ztt" does not fluctuate. 

8.2. Mixed Events Procedure 

The experimental data for iV^/iVs__fluct nations are usually presented in terms of the so 
called dynamical fluctuations 



adyn = sign((T2 - cT^i^.) - ct^,^.| ' , (43) 

where cr^ is defined by Eq. f HTl) . and a^j^ corresponds to the following mixed events 
procedur^. One takes a large number of nucleus-nucleus collision events, and measures 
the numbers of Na and Nb in each event. Then all A and B particles from all events 
are combined into one set. A construction of mixed events is done like the following: 
One fixes a random number = N^ + Ng according to the experimental probability 
distribution P{N), takes randomly particles {A and/or B) from the whole set, fixes 
the values of A^^ and A^^, and returns these A^ particles into the set. This is the mixed 
event number one. Then one constructs event number 2, number 3, etc. 

Note that the number of events is much larger than the number of hadrons, A^, 
in any single event. Therefore, the probabilities pa and ps = 1 — to take the A 
and B species from the whole set can be considered as constant values during the event 
construction. Another consequence of a large number of events is the fact that all A and 
B particles in any constructed mixed event most probably belong to different physical 
events of nucleus-nucleus collisions. Therefore, the correlations between A^^ and Na 
numbers in a physical event are expected to be destroyed in a mixed event. This is the 
main purpose of the mixed events construction. For any function / of Na and A^^ the 
mixed events averaging is then defined as, 



{f{NA,NB))rm. = E ^(^) E /(^^A, A^s) X 
N Na,Nb 

S(N-N,-N,) ^%J^f p1>"B- ■ (44) 

The straightforward calculation of mixed averages can be simplified by introducing 
the generating function Z{x,y), 

Z{x,y) ^ E^iV) E m - Na - Nb)^-^%^^ {xpAf^ {ypBf- 

N Na,Nb 

§ Other dynamical measures, $ |101[|102j and F [SB], can be also used. See the discussion of this point 
in the recent paper |103j . 

II The idealized mixed events procedure appropriate for model analysis is described here. The real 
experimental mixed events procedure is more complicated and includes experimental uncertainties, such 
as particle identification etc. 
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^P{N) {xpA + ypB 



(45) 



TV 



The averages 
One finds: 



are then expressed as x- and y-derivatives of Z{x,y) at x = y = 1. 




{NANB)mix - {Na) 
= paPb ujn {N) . 



\dxdy J 



(46) 
(47) 
(48) 

(49) 



Calculating the Na/Nb fluctuations for mixed events according to (|4T1) one gets: 
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{Na) ' {Nb) ' 

A comparison of the final result in Eq. fl50|) with Eq. fH2|) shows that the mixed event 
procedure gives the same a' for Na/Nb fluctuations as in the GCE formulation for an 
ideal Boltzmann gas, i.e. ua = <^b = ^ and pab = 0. If ojn = 1 (e.g., for the Poisson 
distribution P{N)), one indeed finds o;™^ = w™^ = 1 and Pab = 0. Otherwise, if 
con 7^ 1, the mixed events procedure leads to a;™" ^ 1, o;™^ ^ 1, and to non-zero 
NaNb correlations, as seen from the second line of fl50|) . However, the final result for 
Cmix (Eq]) is still the same simple. It does not depend on the specific form of P{N). 
Non-trivial (c<;™J ^ 1) fluctuations of A^^ and Nb as well as non-zero Pab correlations 
may exist in the mixed events procedure, but they are cancelled out exactly in cr^ix- 



8.3. Hadron Scaled Variances and Correlation Coefficients 

The values of uj.„, ujk, <^p and pxn, PKn, PKn for the HSD simulations of and statistical 
model (SM) in Pb+Pb (Au+Au) central collisions are presented in Figs. and 133 

Let us first comment the SM (see more details in Refs. |98|l99] ). In the SM the scaled 
variances ua and correlation parameter pab approach finite values in the thermodynamic 
limit of large volumes. For central Pb+Pb and Au-|-Au collisions the corresponding 
volumes in the SM are large enough. Finite volume corrections are expected on the level 
of a few percent. We discuss K/ n fluctuations in more details. The n-K correlations pxn 




Figure 38. The SM results in tlie GCE, CE, and MCE ensembles and the HSD resuhs 
(impact parameter b — 0) are presented for the scaled variances w^r, i-^k^ for Pb+Pb 
(Au+Au) collisions at different cm. energies ^/snn- 
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Figure 39. The SM results in the GCE, CE, and MCE ensembles and the HSD 
results (impact parameter 5 = 0) are presented for the correlation parameters pktt, 
PKpi Pp-n for Pb+Pb (Au+Au) collisions at different cm. energies ^/snn 
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are due to resonances having simultaneously K and vr mesons in their decay products. 
In the hadron-resonance gas within the GCE ensemble, these quantum statistics and 
resonance decay effects are responsible for deviations of uk and cj^r from 1, and of pktt 
from 0. The most important effect of an exact charge conservation in the CE ensemble 
is a suppression of the kaon number fluctuation. This happens mainly due to exact 
strangeness conservation and is reflected in smaller CE values of ujk at low collision 
energies in comparison to those from the GCE ensemble. The MCE values of uk and 
are further suppressed in comparison those from the CE ensemble. This is due to exact 
energy conservation. The effect is stronger for pions than for kaons since pions carry a 
larger part of the total energy. An important feature of the MCE is the anticorrelation 
between A^^^ and Nx-, i-e. negative values of pk-k- This is also a consequence of energy 
conservation for each microscopic state of the system in the MCE [63] • The presented 
results demonstrate that global conservation laws are rather important for the values of 
cJtt, ujr) and pk-k- In particular, the exact energy conservation strongly suppresses the 
fluctuations in the pion and kaon numbers and leads to uk < 1 and cJtt < 1 in the MCE 
ensemble, instead of uk > 1 and 0;,^ > 1 in the GCE and CE ensembles. The exact 
energy conservation changes also the n-K correlation into an anticorrelation: instead of 
Pktt > in the GCE and CE ensembles one finds pk-k < in the MCE. The effects of 
global conservation laws and resonance decays are also seen for pxp, Ppn and Up. 

The HSD results for ua and pab shown by the solid lines in Figs. |38] and [39] are close 
to the CE and MCE results for low SPS energies. One may conclude that the infiuence 
of conservation laws is more stringent at low collision energies. The HSD values for 
ua and pab increase, however, at high collision energies and a sizeable deviation of the 
HSD results from those in the MCE SM is observed with increasing energies. 

In Fig. Ho] the HSD scaled variance cu^r for the full acceptance is shown in A+A 
collisions for 6 = 0. The samples of collision events selected experimentally are 3.5% of 
most central collision events in Pb+Pb collisions at the SPS energies and 5% in Au+Au 
collisions at RHIC energies. 

Fig. Ho] presents the HSD results for u^^ in these samples the comparison with the 
HSD results at 6 = 0. One finds much larger values of cj^r in the centrality selected 
samples than for 6 = 0. The effect is especially strong at RHIC energies. Collisions 
with 6 = correspond to up = 0. Thus, in Eq. ( I20l) can be approximately taken as 
c^TT at 6 = 0. The HSD results correspond approximately to up = 0.5 for the 3.5% most 
central Pb+Pb collisions at SPS energies and up = 1 for the 5% most central Au+Au 
collisions at RHIC energies. The results of the MIS f l20|) for are shown by the dashed 
line in Fig. |40]and they are close to the actual values of the HSD simulations for Ut^. 

The fluctuation in the kaon to pion ratio is dominated by the fluctuations of kaons 
alone since the average multiplicity of kaons is about 10 times smaller than that of pions. 
The model calculations of a^Kii) require, in addition to uk, ujn, and pxn values, the 
knowledge of the average multiplicities (Nx) and {N^^). To flx the average multiplicities 
in the SM one needs to choose the system volume. For each collision energy the volume 
of the statistical system has been flxed in a way to obtain the same kaon average 
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Figure 40. (Color online) The HSD results for for Pb+Pb (Au+Au) collisions 
at different cm. energies ^/snn within the full 47r-acceptance. The lower solid line 
corresponds to zero impact parameter {b — 0) and the upper one to the experimentally 
selected samples of collision events. The dashed line reflects the model of independent 
sources (^0)) . 



multiplicity in the SM as in the HSD calculations: {NK)stat = {Nk)hsd- Recall that 
average multiplicities of kaons and pions are the same in all statistical ensembles. The 
SM volume in central Pb+Pb (Au+Au) collisions is large enough and all statistical 
ensembles are thermodynamically equivalent for the average pion and kaon multiplicities, 
since these multiplicities are much larger than 1. 

In Fig. HI] the values of a and amix (in percents) are presented, in the left and right 
panel, respectively, for the SM in different ensembles as well as for the HSD simulations. 
A first conclusion from Fig. UT] (left) is that all results for a in the different models are 
rather similar. One observes a monotonic decrease of a with collision energy. This is 
just because of an increase of the kaon and pion average multiplicities with collision 
energy. The mixed event fluctuations CTmix in the model analysis are fully defined by 
these average multiplicities according to ( l50l) . The values of amix are therefore the 
same in the different statistical ensembles. They are also very close to the HSD values 
because the statistical system volume has been defined to obtain the same kaon average 
multiplicities in the statistical model as in HSD at each collision energy. As seen from 
Fig. HI] (right) the requirement of {NK)stat = {Nk)hsd leads to almost equal values of 
amix in both HSD and the SM. 

Differences between the statistical ensembles as well as between the statistical and 
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Figure 41. Left: The SM results in the GCE, CE, and MCE ensembles as well as 
the HSD results (impact parameter 5 = 0) are presented for cr-100% for K/n ratio in 
Pb+Pb (Au+Au) collisions at different cm. energies ^/snn- Right: The same as in 
the left panel, but for in mixed events defined by (jSOl) 




Figure 42. Left: The results for the K/tt fluctuations at different cm. energies ^/snn 
in the GCE, CE, and MCE ensembles as well as from HSD (impact parameter b = 0) 
are presented for (Tcij,„-100% defined by (|43p . Right: The same as in the left panel but 
forF = aVo-L:.- 



HSD results become visible for other measures of K/tt fluctuations, such as adyn defined 
by fH51) and F = cr^/a^.^. They are shown in Fig. left and right, respectively. 
At small collision energies the CE and MCE results in Fig. |12] demonstrate negative 
values of adyn, respectively F < 1. When the collision energy increases, adyn in the CE 
and MCE ensembles becomes positive, i.e. F > 1. Moreover, the different statistical 
ensembles approach to the same values of adyn and F at high collision energy. In the 
statistical model the values of a and amix approach to zero at high collision energies 
due to an increase of the average multiplicities. The same trivial limit should be also 
valid for adyn in the SM. In contrast, the measure F shows a non-trivial behavior at 
high energies: the SM gives F = 1.05 in high energy limit. The HSD result for F 
demonstrates a monotonic increase with collision energy. An interesting feature of the 
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SM is the approximately equal values of a (and, thus, adyn and F) in the CE and MCE 
ensembles. From Figs. [38] and [39] one observes that both ojk, ^-k and pk-k are rather 
different in the CE and MCE. Thus, as discussed above, an exact energy conservation 
influences the particle scaled variances and correlations. These changes are, however, 
cancelled out in the fluctuations of the kaon to pion ratio. 

It has been mentioned in the literature (see, e.g., |95]) that the fluctuations 
of particle number ratio N^/Nb are independent of volume fluctuations since both 
multiplicities Na and N + B are proportional to the volume. This is not correct. In 
fact, the average multiplicities (Na) and (Nb), but not A^^ and Nb, are proportional 
to the system volume. Let us consider the problem in the SM assuming the presence 
of volume fluctuations with the distribution function Vr(l^) at fixed values of T and 
fj^B- This assumption corresponds approximately to the volume fluctuations in A+A 
collisions because of different impact parameters in each collision event. Under these 
assumptions the SM values of coa, ojb, and pab remain the same for any volume V 
(if only this volume is large enough). However, the average hadron multiplicities are 
proportional to the volume V. Therefore, the SM result for reads, o"^ = ctq Vq/V, where 
Vo is the average system volume, and ctq is calculated with the average multiplicities 
corresponding to this average volume Vq. Expanding Vq/V = Vo/(Vo + SV) in a serious 
of 6V/Vo, one finds to the second order in SV/Vq, 



corresponds to averaging over the volume distribution function Vr(V^) which describes 
the volume fluctuations. As clearly seen from Eq. fl5ip the volume fluctuations influence, 
of course, the K/ vr particle number fluctuations and make them larger. Comparing the 
K/tt particle number fluctuations in, e.g., 1% of most central nucleus-nucleus collisions 
with those in, e.g., 10% one should take into account two effects. First, in the 10% sample 
the average volume Vq is smaller than that in 1% sample and, thus, ctq in Eq. (!5T]) is 
larger. Second, the volume fluctuations (13^ in the 10% sample is larger, and this gives 
an additional contribution to according to Eq. (jSTj). 

One may also consider volume fluctuations at fixed energy and conserved charges 
(see, e.g., Ref. |104] ). In this case the connection between the average multiplicity 
and the volume becomes more complicated. The volume fluctuation within the MCE 
ensemble can strongly affect the fluctuations in the particle number ratios. This 
possibility will be discussed in more detail in a forthcoming study. 

8.4- Excitation Function for the Fluctuations of Ratios 



-^2 _2 1 I 

Cr = CTn 1 + 



(51) 



where 




(52) 



A comparison of the SM results for fluctuations in different ensembles with the data 
looks problematic at present; the same is true for most other theoretical models. This 
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is because of difficulties in implementing the experimental acceptance and centrality 
selection which, however, can be taken into account in the transport approach. In order 
to compare the HSD calculations with the measured data, the experimental cuts are 
applied for the simulated set of HSD events. 

In Fig. Us] the HSD results of adyn for the K/ti, p/tt and K/p ratios are shown 
in comparison with the experimental data by the NA49 Collaboration at the SPS [92] 
and the preliminary data of the STAR Collaboration at RHIC [93l|9ll[in5l|106]. The 
available results of UrQMD calculations (from Refs. [92 t [97 t[T07] ) are also shown by the 
dashed lines. 

For the SPS energies we use the NA49 acceptance tables from Ref. [92]. For 
the RHIC energies we employ the following cuts: in pseudorapidity, |?7| < 1, and in 
transverse momentum, 0.2 < pt < 0.6 GeV/c for kaons and pions and 0.4 < pt < 
1 GeV/c for protons [93 t [M tfT05lll06] . We note also, that the HSD results presented in 
Fig. Us] correspond to the same centrality selection as in the experiment: the NA49 data 
correspond to the 3.5% most central collisions selected via the veto calorimeter, whereas 
in the STAR experiment the 5% most central events with the highest multiplicities in 
the pseudorapidity range |?7| < 0.5 have been selected. 

One sees that the UrQMD model gives practically a constant ci^^, which is by 
about 40% smaller than the results from HSD at the lowest SPS energy. This difference 
between the two transport models may be attributed to different realizations of the 
string and resonance dynamics in HSD and UrQMD: in UrQMD the strings decay first 
to heavy baryonic and mesonic resonances which only later on decay to 'light' hadrons 
such as kaons and pions. In HSD the strings dominantly decay directly to 'light' hadrons 
(from the pseudoscalar meson octet) or the vector mesons p, u and K* (or the baryon 
octet and decouplet in case of baryon number ±1). Such a 'non-equilibrated' string 
dynamics leads to stronger fluctuations of the K/n ratio. Note that all differences 
between SM and transport models, as well as between different versions of the transport 
models, become clearly seen at the lower bombarding energies. This is only because 
of using adyn as a measure of the K/ti ratio fluctuations. If one uses F = cr^/cj^-^ as 
a measure of the K/n fluctuations the conclusion will be opposite: as Fig. |12] (right) 
demonstrates the difference between the SM and HSD predictions measured in F would 
increase with collision energy. 

The excitation function of for the K/p ratio is presented in the middle panel of 
Fig.HHl The 'stars' correspond to the preliminary STAR data in Au+Au collisions [106j . 
The HSD results for show a weak energy dependence in both SPS and RHIC energy 
regions. However, as follows from our analysis this observable is very sensitive to the 
experimental acceptance. In order to illustrate it, we extended the HSD calculations 
for Au+Au within the STAR acceptance for SPS energies also - cf. the upper red line 
with solid dots, which it is higher then the HSD result for Pb+Pb obtained for the 
NA49 acceptance - lower blue line with open dots. The influence of the experimental 
acceptance is clearly seen at 158 AGeV where a switch from the NA49 to the STAR 
acceptance leads to a jump in by more than 3%. On the other hand, our calculations 
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Figure 43. The HSD results for the excitation function in adyn for the K/n, 
K/p, p/tt ratios within the experimental acceptance (solid line) in comparison to the 
experimental data measured by the NA49 Collaboration at SPS WI} and by the STAR 
Collaboration at RHIC [93ll94|jl05lll06) . The UrQMD calculations are shown by dotted 
lines. 
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for Pb+Pb (3.5% central) and for Au+Au (5% central) collisions - performed within the 
NA49 acceptance for both cases at 160 AGeV - shows a very week sensitivity of on 
the actual choice of the collision system and centrality - cf. the coincident open circle 
and triangle at 158 AGeV in the middle panel of Fig. HS] Indeed, new data from NA49 
and STAR should clarify the situation. 

At SPS energies the HSD simulations lead to negative values of adyn for the proton 
to pion ratio - cf. Fig. HSl lower panel. This is in agreement with the NA49 data 
in Pb+Pb collisions. On the other hand HSD gives large positive values of at 
RHIC energies which strongly overestimate the preliminary STAR data for Au+Au 
collisions |105] . 

9. Summary 

The present review has been devoted to a systematic study of fluctuations and 
correlations in heavy-ion collisions within the HSD and UrQMD transport models. This 
provides a powerful tool to simulate experimental collisions on an event-by-event basis. 

In Section [3] the centrality dependence of multiplicity fluctuations has been studied 
in fixed target A+A experiments. The centrality selection in Pb+Pb collisions at 
158 AGeV has been performed in full correspondence to the NA49 experimental 
procedure: only the samples of events with fixed numbers of the projectile participants, 
N'p°\ have been considered. A decrease of the scaled variance up for the participant 
number fluctuations with Np^"-' in central A+A collisions has been found. The transport 
model results for the scaled variances Ui of negative, positive, and all charged hadrons 
in Pb+Pb at 158 AGeV have been presented and interpreted. The samples with 
jyproj _ 20 + 30 show a maximum of fluctuations of the number of target nucleon 
participants Np^^. The final hadron multiplicity fluctuations exhibit an analogous 
behavior. This explains the large values of the HSD and UrQMD scaled variances Ui for 
positive, negative, and all charged hadrons in the target hemispheres and in the full Air 
acceptance. The asymmetry between the projectile and target participants - introduced 
in the data samples by the trigger condition of fixed Np'"-' - can be used to explore 
different dynamical scenarios for A+A collisions by measuring the final multiplicity 
fluctuations as a function of rapidity (cf. Fig. |9]). The analysis reveals that the NA49 
data indicate a rather strong mixing of the longitudinal flows of the projectile and 
target hadron production sources. This sheds new light on the A+A reaction dynamics 
at top SPS energies. The data demonstrate a significantly larger amount of mixing 
than generated in hadron-string transport approaches. This probably points towards 
the presence of strongly interacting quark-gluon phase. 

In Section H] the particle number fluctuations in central A+A collisions from 10 to 
21300 AGeV have been studied within the HSD transport model. HSD predicts that the 
scaled variances Ui in central A+A collisions remain close to the corresponding values in 
p + p collisions and increase with collision energy with the multiplicity per participating 
nucleon, Ui oc rii. The scaled variances Ui calculated within the statistical hadron- 
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resonance gas (HG) model along the chemical freeze-out line show a rather different 
behavior: the Ui approach finite values at high collision energy. At the top RHIC 
energy ^/sNN = 200 GeV the HSD results for a;i(HSD) are already about 10 times larger 
than the corresponding micro-canonical (MCE) HG values for a;j(MCE). Thus, the HSD 
and HG scaled variances Ui show a different energy dependence and are very different 
numerically at high energies. However, a comparison with the NA49 data of very central, 
< 1%, Pb+Pb collisions in the SPS energy range does not distinguish between the HSD 
and MCE HG results. This occurs because of two reasons: First, the MCE HG and 
HSD results for Ui at SPS energies are not too much different from each other at SPS 
energies as well as from the Ui values in p-|-p collisions. Second, small experimental 
values of the acceptance, q = 0.04 -^ 0.16, and about 10% possible ambiguities coming 
from the acceptance scaling relation ( l28l) make the difference between the HSD and 
MCE HG results almost invisible. New measurements of Ui for samples of very central 
A-l-A collisions with large acceptance at both SPS and RHIC energies are needed to 
allow for a proper determination of the underlying dynamics. 

In Section |5] the event-by-event multiplicity fluctuations in A+A collisions have 
been studied for different energies and system sizes within the transport approach. 
This study is in full correspondence to the future experimental program of the NA61 
Collaboration at the SPS. The central C-l-C, S+S, In+In, and Pb-|-Pb nuclear collisions 
from Eiab= 10, 20, 30, 40, 80, 158 AGeV have been investigated. The influence of 
participant number fluctuations on hadron multiplicity fluctuations has been emphasized 
and studied in detail. To make these 'trivial' fluctuations smaller, one has to consider the 
most central collisions. Indeed, one needs to make a very rigid selection - 1% or smaller 
- of the 'most central' collision events. In addition, one wants to compare the event-by- 
event fluctuations in these 'most central' collisions for heavy and for light nuclei. Under 
these new requirements different centrality selections are not equivalent to each other. 
As a consequence, there is no universal geometrical selection by the impact parameter. 
This is a new and serious problem for theoretical models (e.g. for hydrodynamical 
models) in a precision description of the event-by-event fluctuation data. The above 
statements have been illustrated by the 6 = selection criterium. For light nuclei 
even these 'absolutely central' geometrical collisions lead to rather large fluctuations 
of the number of participants, essentially larger than in the 1% most central collisions 
selected by the largest values of the projectile participants Np'"-' . The multiplicity 
fluctuations calculated in these samples show a much weaker dependence on the atomic 
mass number A than for the criterium of zero impact parameter 6 = 0. A monotonic 
energy dependence for the multiplicity fluctuations is obtained in both the HSD and 
UrQMD transport models. The two models demonstrate a similar qualitative behavior 
of the particle number fluctuations. Our study has also demonstrated a sensitivity of 
the multiplicity fluctuations to some specific details of the transport models. Thus, the 
present HSD and UrQMD results for the scaled variances provide a general trend of their 
dependencies on A and Eiab and indicate quantitatively the systematic uncertainties. It 
has to be stressed again, that HSD and UrQMD - in their present formulations - do not 
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include explicitly a phase transition to the QGP. The expected enhanced fluctuations - 
attributed to the critical point and phase transition - can be observed experimentally on 
top of a monotonic and smooth 'hadronic background'. The most promising signature 
of the QCD critical point would be an observation of a non-monotonic dependence of 
the scaled variances with bombarding energy Eiab for central A+A collisions with flxed 
atomic mass number. In the flxed target SPS experiments the measurements of a;_, 
a;+, and Uch are preferable in the forward hemispheres. In this case the remaining 
small fluctuations of the number of target participants N'^^^ in the 1% most central 
collisions - deflned by the number of the projectile participants Np'"-' - become even 
less important, as they contribute mainly to the particle fluctuations in the backward 
hemisphere. These flndings are expected to be helpful for the optimal choice of collision 
systems and collision energies for the experimental search of the QCD critical point. 

In Section [6] the HSD transport model has been used to calculate the scaled 
variance of participant number fluctuations, up, and the number of z'th hadrons per 
nucleon accepted by the mid-rapidity PHENIX detector, qiUi, in different Beam-Beam 
Counter centrality classes. The HSD model for N+N collisions has been also used to 
estimate the fluctuations from a single source, u*. It has been found that the model 
of independent sources based on these HSD results leads to a good agreement with the 
PHENIX data [761177]. In different 5% centrality classes, up goes down and g^nj goes up 
with increasing of (Np). This results in a non-monotonic dependence of u'^^'^ on (Np) as 
seen in the PHENIX data. Thus, one may conclude that the centrality dependence of 
the fluctuations seen in the present PHENIX data are the consequences of participant 
number fluctuations. To avoid a dominance of the participant number fluctuations we 
stress again that one needs to analyze most central collisions with a much more rigid 
(< 1%) centrality selection. 

The system size event-by-event fluctuations induce also rapidity forward-backward 
correlations in relativistic A+A collisions. The analysis has been based on two 
independent models - the microscopic HSD transport approach and a 'toy' wounded 
nucleon model realized as a Glauber Monte Carlo event generator. It has been shown 
that strong forward-backward correlations arise due to an averaging over many different 
events that belong to one 10% centrality bin. In contrast to average multiplicities, the 
resulting fluctuations and correlations depend strongly on the speciflc centrality trigger. 
For example, the centrality selection via impact parameter b (used in most theoretical 
calculations) and via N^j^'^ (used experimentally) lead to rather different values of 
participant number fluctuations up and forward-backward correlation parameters p. 
These different centrality selections give also quite different dependencies of up and p 
on centrality. The HSD simulations reveal strong forward-backward correlations and 
reproduce the main qualitative features of the STAR data in Au+Au collisions at RHIC 
energies [HSIIHS]. The forward-backward correlations can be studied experimentally 
for smaller size centrality bins deflned by N^^-^ . When the size of the bins decreases, 
the contribution of 'geometrical' fluctuations should lead to weaker forward-backward 
correlations and to a less pronounced centrality dependence. Note, that the 'geometrical' 
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fluctuations discussed here are in fact present in all dynamical models of A+A collisions. 
Thus, they should be carefully accounted before any discussion of new physical effects 
is addressed. A future experimental analysis - in the direction examined here - should 
clarify whether the observed correlations by the STAR Collaboration at RHIC contain 
really additional contributions from 'new physics'. 

In Section [7]the transport approach has been used to investigate the event-by-event 
fluctuations of baryon number and electric charge in A+A collisions. The study has been 
based on the microscopic HSD transport model which allows to analyze on the event- 
by-event basis the influence of two main sources of the conserved charge fluctuations - 
the fluctuations of the number of participants and the finite experimental acceptance. 
It has been found that the fluctuations in the number of target participants strongly 
influences the baryon number and electric charge fluctuations. The consequences 
depend crucially on the dynamics of the initial flows of the conserved charges. For 
a better understanding of the HSD results three limiting groups of models for A+A 
collisions have been discussed: transparency, mixing and reflection. These "pedagogical" 
considerations indicate that the HSD model shows only a small mixing on initial baryon 
flow and is closer to the transparency models. This makes important the findings of 
Refs. [26 | lis t [T08| 1109] about the influence of the partonic degrees of freedom on the 
initial phase dynamics. This phase might increase the mixing of conserved charges by 
additional strong parton-parton interactions. 

In Section |8] the event-by-event multiplicity fluctuations of pions, kaons, protons as 
well as their correlations and ratio fluctuations in central Au+Au (or Pb+Pb) collisions 
from low SPS up to top RHIC energies have been studied within the HSD transport 
approach and in the statistical hadron-resonance gas model for different statistical 
ensembles - the grand canonical ensemble (GCE), canonical ensemble (CE), and micro- 
canonical ensemble (MCE). The HSD results at SPS energies are close to those in 
the CE and MCE statistical models. This indicates a dominant role of resonance 
decays and global conservation laws for low energy A+A collisions. On the other 
hand, substantial differences in HSD and statistical model results have been observed at 
RHIC energies which can be attributed to non-equilibrium dynamical effects in the HSD 
simulations. Thus, the scaled variances of different particle species may serve as good 
observables to probe the amount of equilibration achieved in central A+A collisions at 
RHIC energies. On the other hand, it has been found that the observable (Jdyn, which 
characterizes ratio fluctuations, appears to be rather sensitive to the details of the 
model at low collision energies. It has been found that HSD can qualitatively reproduce 
the measured excitation function for the K/tt ratio fluctuations in central Au+Au (or 
Pb+Pb) collisions from low SPS up to top RHIC energies. The HSD results for cr^^^ 
appear to be close to the NA49 data at the SPS. The data for a^J^ in Pb+Pb collisions 
at the SPS energies will be available soon and allow for further insight. A comparison of 
the HSD results with preliminary STAR data in Au+Au collisions at RHIC energies are 
not fully conclusive: adyn from the HSD calculations is approximately in agreement with 
data |106] for the kaon to proton ratio, but overestimates the experimental results |105] 
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for the proton to pion ratio. New data on event-by-event fluctuations in Au+Au at 
RHIC energies will help to clarify the situation. 

It has been argued in the Introduction that event-by-event fluctuations in nucleus- 
nucleus collisions can be considered as a probe for the phase transition and the QCD 
critical point. The transport models have been employed to study the fluctuations and 
correlations including the influence of experimental acceptance as well as the centrality 
selection, the dependence on the system size and collision energy. 

So far, the transport models employed do not include an explicit phase transition 
to the QGP as well as a dynamical description of hadronization. Recent steps have 
been taken in this direction within Parton-Hadron-String-Dynamics (PHSD) approach 
|110yillj that incorporates the partonic dynamics in line with an equation of state from 
lattice QCD. Nevertheless, the comparison of results from microscopic transport models 
with experimental data will allow to separate the effects indicated by a phase transition 
to the QGP and the QCD critical point. 
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